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Preface 


This  volume,  the  second  of  three  volumes  which  examine  the  noise 
attenuation  value  of  a  flexibly-supported  bearing,  is  divided 
into  two  parts! 

Part  I  -  Attenutatlon  of  Rotor  Unbalance  Forces 
jBy  Flexible  Bearing  Supports 

Part  XI  -  Unbalance  Response  Of  A  Uniform  Elastic 

Rotor,  Supported  In  Ibamped  Flexible  Bearings 

The  other  two  volumes  complete  a  study  aimed  specifically  at  an 
investigation  into  the  effect  of  a  hydraulically  supported, 
tilting-pad.  Journal  bearing  on  the  attenuation  of  noise  origi¬ 
nating  from  rotor  unbalance.  Volume  1,  Spring  and  Damping 
Coefficients  Por  The  Tllting-Pad  Journal  Bearing,  provides  an 
analytical  method  7 or  determining  the  spring  ana  damping 
properties  of  the  bearing  oil-film  with  the  results  presented 
iri  curves  for  typical  tilting-pad  bearing  geometries.  The 
third  volume  includes  (1)  a  generalized-rotor  analysis  and 
(2)  experimental  results. 

Mechanical  Technology  Incorporated  waB  primarily  responsible 
for  the  analytical  portion  of  this  Btudy  -  while  Westlnghouse 
Electric  Corporation  designed  and  conducted  the  experimental 
test. 
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ABSTRACT 


This  report  is  concerned  with  the  attenuation  of  the  force  transmitted 
by  an  unbalanced  rotor.  The  attenuation  is  achieved  by  a  flexible  bearing 
support.  The  report  presents  an  analysis  of  the  unbalance  vibrations 
of  a  flexible,  symmetrical ,  two-maso  rotor  supported  in  fluid  film  bearings 
which  in  turn  are  mounted  in  flexible  supports.  The  fluid  film  in  the 
bearings  possesses  both  flexibility  and  damping.  The  analysis  takes  into 
account  both  static  and  dynamic  rotor  unbalance  and  gives  results  for  the 
rotor  amplitudes  and  the  transmitted  force  as  functions  of  the  system 
parameters.  The  analysis  has  been  programmed  for  a  digital  computer  and 
a  description  of  the  two  computer  programs  is  included.  Numerical  results 
have  been  obtained  and  are  summarized  in  24  design  charts. 


INTRODUCTION 


Among  the  principal  noise  sources  on  board  a  ship  are  the  vibrations  trans¬ 
mitted  from  the  various  pieces  of  rotating  machinery  such  as  the  main  pro¬ 
pulsion  unit,  the  auxiliary  machinery,  etc.  In  attempting  to  reduce  the 
transmitted  vibrations  two  general  approaches  are  available: 

a)  eliminate  the  causes  responsible  for  the  generated  noise  by 
better  rotor  balancing,  closer  manufacturing  tolerances,  etc. 

b)  attenuate  the  noise  by  means  of  vibration  isolation 

Experience  has  shown  that  the  former  approach  yields,  at  best,  moderate 
gains  due  to  the  limitations  imposed  by  practical  considerations.  Hence, 
increased  attention  is  given  instead  to  the  possible  methods  of  vibration 
isolation.  A  key  consideration  in  the  application  of  this  approach  is  to 
attenuate  the  noise  as  close  to  its  source  as  feasible  and  prevent  the  vibra¬ 
tions  from  setting  too  heavy  masses  in  motion.  This  may  be  accomplished  by 
mounting  the  bearing  housings  on  flexible  supports. 

It  is  the  purpose  of  this  report  to  present  an  analysis  of  the  noise  atten¬ 
uating  characteristic  of  a  flexible  bearing  support.  Specifically  the  analysis 
establishes  the  equations  for  determining  the  support  stiffness  which  will 
achieve  a  desired  noise  attenuation.  The  analysis  takes  into  account  the 
flexibility  of  the  rotor,  the  damping  and  stiffness  characteristics  of  the 
bearing  film,  the  mass  of  the  bearing  housing  and  the  stiffness  and  damping 
of  the  support.  Both  static  and  dynamic  rotor  unbalance  is  considered. 

Numerical  results  have  been  computed  for  a  wide  range  of  support  and  rotor 
parameters.  The  results  are  presented  in  24  graphs  giving  the  transmitted 
force  and  the  rotor  amplitude  as  functions  of  the  rotor  speed.  The  graphs  are 
intended  to  be  used  for  selecting  the  support  parameters  (e.g.  stiffness)  to 
achieve  a  desired  noise  attenuation  over  a  given  speed  range. 
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DISCUSSION 

The  unbalance  which  is  always  present  in  a  rotor  gives  rise  to  a  dyna¬ 
mical  force  at  the  bearings  with  the  same  frequency  as  the  speed  of  the 
rotor.  Some  attenuation  of  the  force  takes  place  in  the  bearings  since 
the  fluid  film  possesses  both  flexibility  and  damping  (see  Refs  1,2,3,  4 
and  5).  However,  to  obtain  an  appreciable  force  attenuation  it  becomes 
necessary  to  vibration  isolate  the  bearings  by  means  of  a  flexible  support 
(e.g.  a  hydraulic  support).  The  magnitude  of  the  resulting  attenuation  depends 
on  the  stiffness  of  the  support  and  the  dynamical  characteristics  of  the 
rotor-bearing  system.  It  is  the  purpose  of  the  present  report  to  study 
this  attenuation  and  to  determine  the  effect  of  the  parameters  of  the 
system 

The  system  is  represented  by  the  model  in  Fig.  1.  It  consists  of  a  flexible, 
symmetrical,  two-mass  rotor  supported  in  two  journal  bearings.  Having  two 
rotor  masses  makes  it  possible  to  consider  both  static  and  dynamic  rotor 
unbalance.  Each  bearing  is  mounted  on  a  flexible  support  with  a  specified 
stiffness  and  a  specified  damping.  The  journal  bearing  is  characterized  by 
4  spring  coefficients  and  4  damping  coefficients  derived  from  lubrication 
theory  (Ref s . 1 , 2 , 3 ,4 , 5) .  The  mass  of  the  bearing  housing  is  also  included 
since  it  affects  the  force  transmission. 

On  the  basis  of  the  selected  model  the  equations  of  motion  are  set  up  for  the 
rotor  and  the  support.  Because  of  the  large  number  of  parameters  it  is  of 
limited  value  to  derive  a  closed  form  solution.  Instead,  the  equations  are 
reduced  to  a  form  convenient  for  numerical  evaluation  and  programmed  for  a 
digital  computer.  Two  computer  programs,  both  for  the  IBM  1620  computer,  have 
been  written.  They  are  described  in  detail  in  the  Appendix  including  instruc¬ 
tions  for  using  the  programs.  The  results  from  the  programs  include  the  force 
transmitted  to  the  foundation,  the  rotor  amplitudes  and  the  amplitude  of  the 
bearing  housing.  Extensive  calculations  have  been  performed  and  the  data  have 
been  plotted  in  Figs. 2  to  25.  The  employed  values  of  rotor  stiffness  and  weight 
and  of  support  stiffness  are  selected  to  cover  the  range  normally  encountered  in 
Navy  applications.  The  use  of  the  graphs  is  explained  in  the  following  section. 


A  study  of  the  graphs  reveals  that  a  very  significant  attenuation  of  the 
transmitted  force  can  be  obtained,  at  least  in  theory  Note  that  the 
"unattenuated"  force  is  not  shown  in  the  graphs  but  would  appear  as  a  series 
of  straight  lines.  They  are  determined  by  a  simple  relationship  given  in  the 
following  section.  However,  even  small  amounts  of  damping  in  the  support 
diminish  the  attenuation  and  in  practice  this  is  unavoidable.  Furthermore, 
system  resonances  tend  to  become  accentuated  in  a  way  not  found  in  the  con¬ 
ventional  construction,  e  g.  the  resonance  of  the  bearing  housing  becomes 
important.  Therefore,  the  graphs  should  serve  as  a  guide  line  only,  they  are 
not  intended  to  be  final  design  charts. 

Since  it  is  known  that  a  flexible  support  may  adversely  affect  the  stability 
of  the  bearing,  i.e.  the  speed  ar  onset  of  oil  whip  may  be  lowered,  a  study 
of  the  stability  is  also  undertaken.  It  is  concluded  that  even  if  a  flex¬ 
ible  support  does  lower  the  oil  whip  speed  when  there  is  no  damping  in 
the  support  it  takes  only  a  rather  small  amount  of  damping  to  restore,  and 
even  increase  the  stability  limit. 


RESULTS 


The  main  purpose  of  the  numerical  results  (Fig. 2  to  25)  is  to  illustrate 
the  attenuation  of  the  transmitted  force  due  to  a  flexible  support.  To 
obtain  a  high  attenuation  as  required  in  Naval  applications  the  support 
is  very  soft  in  comparison  with  the  bearing  stiffness. 

The  rotor  configuration  is  shown  in  Fig.  1.  The  rotor  is  symmetric  and 
consists  of  two  masses  M  on  a  flexible  shaft.  There  are  two  bearings  which 
are  flexibly  supported  with  the  support  stiffness,  Kp .  In  addition,  the  mass 
m  of  the  bearing  housing  is  included  because  the  support  resonance  may  be¬ 
come  important  when  the  stiffness  is  small.  The  support  damping  has  been 
neglected  since  it  is  in  general  kept  small  in  order  to  achieve  the  force 
attenuation. 


The  results  are  given  in  form  of  dimensionless  parameters: 


dimensionless  transmitted  force: 


dimensionless  journal  amplitude: 


tci 

Wtf 

U * 


Xj 

dimensionless  rctor  mass  amplitude:  Tr 
speed  ratio: 


CO 


The  rotor-bearing-support  parameters  are; 


CK»x 

W 

CcuCxx 

w 

C  Kp 

W 


bearing  stiffness: 

bearing  damping: 

support  stiffness: 

.  hi 

support  mass:  t^pj1 

Wot 

rotor  flexibility  parameter;  0  =  _ 
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where: 


w 

M 

ck 

r 

6 

C 

U) 

hi 

KP 
K 
C 
r 

X. 
X, 


XX 


XX 


bearing  reaction,  lbs. 

2 

half  of  the  total  rotor  mass,  lbs. sec  /in 

influence  coefficient  for  rotor  (inverse  rotor  stiffness), 
Equation  (4) ,  in 
lbs. 

distance  between  rotor  masses  divided  by  rotor  length, 

Equation  (5). 

eccentricity  of  rotor  masses,  inch 
radial  bearing  clearance,  inch 
rotor  speed,  rad/sec. 

I/VMoI  J  critical  speed  of  rigidly  supported  rotor,  rad/sec. 

2 

support  mass,  lbs. sec  /in 

support  stiffness,  lbs/in 

bearing  stiffness,  lbs/in 

bearing  damping  coefficient,  lbs.sec/in 

force  transmitted  to  foundation,  lbs. 

amplitude  of  rotor  mass,  inch 

amplitude  of  journal,  inch. 


The  numerical  data  cover  the  following  range  of  rotor-support  parameters: 

?=  cf*  =  .1 ,  .3,  1,3  / 0 

-  I0"3,  3  *  jo-3  ,  ICf2  3  *  )0'2 

Since  in  the  present  case  the  bearing  is  much  stiffer  than  the  support) it  has 
virtually  no  influence  on  the  results.  For  completeness  the  bearing  is  assigned 
the  coefficients: 

=  3.75  «  4.1  o 

which  is  representative  of  a  4-shoe  tilting  pad  journal  bearing. 


Interpretation  of  the  Charts 


The  data  apply  to  both  static  and  dynamic  unbalance  such  that  for  static 
unbalance  £  -/  and  o(  corresponds  to  the  first  flexural  rotor  critical  speed, 
whereas  for  dynamic  unbalance  f  represents  the  distance  between  rotor  masses 
and  d  corresponds  to  the  second  flexural  rotor  critical  speed. 


In  Figs.  2  to  25  two  resonances  are  evident.  The  first  resonance  is  the 
rotor  resonance  which  for  a  soft  support  occurs  at: 


(—) 

twilit*  rti»hA*tt 

or  in  dimensional  form: 


LO  . 

tyttw,  reio^Anct 


taAiiai 

s  tc 


rttni'L 

ftc 


C  1st  ) 


(Zhd  IVipi If) 


Note  that  for  the  second  mode  the  angular  restoring  stiffness  is  kl\ 

{  t  -  span  between  bearings)  and  the  transverse  mass  moment  of  inertia  of  the 
rotor  is  leading  to  the  stated  result  for  the  resonant  speed. 


The  second  resonance  shown  by  the  graphs  is  the  resonance  of  the  support: 

an 


ISnffrt  rtSththt, 

or  in  dimensional  form: 


COr 


—  ■/ ItKp  ^  »/  — — 

~  Yo (hi  fee 


Skfpcrt  irtiokihtt 

where  represents  the  rotor  stiffness.  Hence,  in  schematic  form  Figs.  2 


to  25  become: 
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TrtinshuHeJ 
font : 

UE 

w<$ 


‘3our*wJ  At»ip li'-tucfe  : 
or 

X 

floUr  ln«<S  ^"l 


Sfetd  -  ^ 


Sp*"* ;  &>„ 


Since  the  support  dan  ping  has  been  neglected  the  resonance  values  are  very 
large  in  the  graphs.  However,  some  damping  is  always  present  in  the  support 
and  an  estimate  of  the  peak  amplitude  and  transmitted  force  at  the  rotor 


resonance  can  be  obtained  from: 


peak  journal  amplitude: 


peak  transmitted  force: 


(l£f) 

\  W<f  ip 


CKe  fix,) 

W  '  <f  'f**k 


In  dimensional  form: 


journal  amplitude: 

incK 

d 

P  = 

transmitted  force:  r  p«*k 

Kp 

Us 

where: 

d  -  support  damping,  lbs.sec/ir. 


When  the  rotor  is  considered  rigid  and  the  bearings  are  also  taken  to 
be  rigid  the  transmitted  force  becomes: 


F  -  f  lbs.  (rigid  rotor) 


—  1  ( unuliehuaUA  font) 

Comparing  this  value  with  the  values  of  the  transmitted  force  for  a 
flexible  support  the  magnitude  of  the  force  attenuation  can  be  estimated. 
The  above  "unatt.enuated"  force  would  appear  as  straight  lines  with  a 
slope  of  2  in  Figs.  2  to  9.  There  would  be  a  line  for  each  rotor  flexi¬ 
bility  parameter. 


or  in  dimensionless  form: 

icf 
w  i 
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Numerlcal  Example 

Let  it  be  desired  to  vibration  isolate  a  steam  turbine  rotor  with  the 
following  data: 

Total  weight:  4,400  lbs. 

6  2 

Transverse  mass  moment  of  inertia  around  CG:  I  =  1.35-10  lbs. in 
Span  between  bearings:  =  74  inch. 

First  flexural  critical  speed:  =  25,200  RPM 

Second  flexural  critical  speed:  =  73,000  RPM 

The  rotor  is  approximately  symmetrical  around  the  CG  so  that: 

Bearing  reaction:  W  =  2,200  lbs 

2 

Rotor  mass  per  bearing:  M  =  2,200  lbs  =  5.7  lbs. sec  /in. 


From  the  formula: 


the  influence  coefficient  et  is  calculated  as: 


ID* 
ir*  i* 


2.5 ‘(0  (1st  mode) 


[  3*10  1  (2nd  mode) 


For  the  second  mode  the  two  rotor  masses  M  are  separated  by  the  distance  ( ) 
where  l  is  the  rotor  span.  Thus: 

M  -  I 

from  which:  . 

*•  -  J_\/  =  1 75- 

5  W  Z2O0 

For  the  first  mode  £  =  1  by  definition. 

Setting  the  radial  bearing  clearance: 

-3 

C  =  4.10  inch 


-li¬ 


the  rotor  flexibility  parameter  ^  becomes: 

h'il’IO  *  (1st  mode) 
7-3*/0  (2nd  mode) 


These  values  are  outside  the  range  used  in  the  graphs  making  it  necessary 
to  employ  a  scale  factor  on  the  influence  coefficient.  This  can  be  done 
because  the  rotor  is  very  stiff.  Hence,  if  e<  is  multiplied  by  7.3  and 
13.7  respectively  the  new  values  becomes  J  =.l  for  both  the  first  and  second 
mode.  The  corresponding  flexural  critical  speeds  become: 

25,200 

£U,  =  =  9,300  RPM  (by  scaling) 

~  y '{ 3  ~ 7  =  *-9,700  RPM  (by  scaling) 

In  order  to  select  the  support  stiffness  Kp  it  is  seen  from  Figs.  2  to  9 
that  it  is  necessary  to  specify  the  lowest  speed  at  which  attenuation  is  re¬ 
quired.  Then  Kp  must  be  selected  such  that  the  system  resonances  are  well 
below  this  specified  speed.  In  the  present  case  it  is  desired  to  achieve  a 
substantial  attenuation  at  1100  RPM  or  for  use  in  the  charts: 


Id  _  i^ioo 

Co,  9,300 

CJ  _  1,100* 

0)z  ~  19,700 


.118 

.056 


The  "unattenuated"  transmitted  force  becomes  a  straight  line  in  Figs.  2  to  9 
determined  by: 


I8L,  ■  *(Sf  ■  »(S) 


/u*Att«hfcAtc<(  ? 

i.e.  a  line  passing  through  two  points: 


(lo']lo'%)  (l  0~\l0~‘) 
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It  is  found  that  =  10  ^  is  the  highest  value  of  the  support  stiffness 

which  gives  attenuation  for  both  modes  at  1,100  RPM.  From  Fig.  6  (or  7); 


actual 

(fii  \  force 

M  rcF 

for  l.lOOrpm  W<f 
1st  mode  (static  unbalance)  .118  .01 

2nd  mode  (dynamic  unbalance)  .056  .0135 


"unattenuated"  Attentuation 


(&), 


force 


w tf/estliti.. 

.135 

.03 


at  1,100  RPM 


22.6 

6.9 


The  actual  value  of  the  support  stiffness  IC, 

.  .  a  a  .  .  r 


Kp  -  ^  lo"*  = 


2,200 
tl^l° 


The  bearing  housing  weighs  100  lbs.  Thus: 


is  calculated  as: 

& 


W)  100 

f2M  = 


(1st  mode) 
(2nd  mode) 


Hence,  Fig.  6  applies  for  the  first  mode  and  Fig.  7  for  the  second  mode. 

However,  for  the  original  calculated  values  of  the  rotor  flexibility  parameter 
the  support  resonance  is  very  insignificant. 

The  journal  and  rotor  amplitudes  can  be  found  from  Figs.  14  and  22  (or  Figs.  15 
and  23  for  the  second  mode).  Since  the  rotor  is  very  stiff  there  is  no  difference 
between  the  two  amplitudes  in  the  speed  range  of  interest,  say  up  to  10,000  RPM. 
Still  using  the  scale  factor  on  o(  the  rotor  resonances  occur  at: 


f  215  PPM  (1st  mode) 
\  630  RPM  (2nd  mode) 


Let  the  rotor  unbalance  be  4  oz.inch  per  rotor  mass,  i.e.  the  rotor  eccentricity 
6  becomes: 

At  speeds  above  1200  to  1300  RPM  the  journal  amplitude  is  practically  constant: 


,H3'/0  *  (1st  mode) 
k  ,24’<0  I *»d i  (2nd  mode) 


Let  there  by  a  small  damping  in  the  support,  given  by  the  damping  coefficient  & 
and  assume  the  value: 

d-  20 


Hence,  the  peak  amplitude  at  resonance  can  be  calculated  from  the  previously 
given  formula: 


5525olf TvoT~ 

2200  60  ll Woo 


MS  (1st  mode) 
^•2  (2nd  mode) 


In  actual  values  the  peak  journal  amplitude  becomes: 


ANALYSIS 


In  order  to  account  for  both  the  first  and  the  second  critical  speed 
and  also  consider  both  static  and  dynamic  unbalance  the  rotor  is  re¬ 
presented  by  a  symmetrical  two-mass  rotor: 


M  M 


Rotor  Equations 

Let  the  rotor  stiffness  be  described  by  the  influence  coefficients  ol«« 
and  ,  applying  at  the  concentrated  masses.  The  rotor  unbalance 

is  introduced  by  giving  each  mass  an  eccentricity  6  .  Hence,  the 
equations  of  motion  are 

(i)  X,  =  ("  Mx,  +  McoVcoSfot)  +  otoi,^  My,  t  McoVreswt)  + 


where  the  upper  sign  applies  to  static  unbalance  (first  mode)  and  the 
lower  sign  is  for  the  dynamic  unbalance  (second  mode).  Eq.(l)  may  also 
be  written: 


(2)  =-oCMx,+ol  McoVfoSCot 

(3)  51  ~  f  f2  *  ~  o(  +  at  Me oV  Sihtot 


with  the  convention: 


(4) 

(5) 


oL  = 


J 


static  Uhb*lahf</  (irfi  friotfc 
dtjnar*> 'c  uhbAl«hf^  SecohJ 
Static  U«fcal«iicc(  ■fii-st  modi 
dijnavtM'c  Uhb*(<?h ctt  Stcet\A 
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The  rigid  support  critical  speed  is  then  given  by: 


(6) 

4>r:Vi 

Set: 

(7) 

x,  =  (Pe  { f  e!~l 

(8) 

X,  =  «  t  ¥  e'i 

(9) 

*3=  e;“‘f 

and  similarly 

for  the  y-direction. 

to  get: 

(10) 

./  _  *2  +  (wj 

i  -  rar 

(ID 

U  =  -  i  (&,)* 

* 

Bearing  Equations 

The  rotor  is  supported  in  two  identical  bearings.  Each  bearing  is  re¬ 
presented  by  four  spring  coefficients:  K*x f  ;  and  four 

damping  coefficients:  )  Cxtj  ;  C^x  f  ^  ,  A  force  balance  yields: 

(12)  ^  MX|  +  MbJ  dcoSU)t\ ~  £  fo-f*]-  Ktx  ^-Xj)  +  CyxfXj-XjJ'f"  kx^  ( 

(13)  ^[-My,+Mu>Vsma]=  |  + 


Set: 

(14) 


\s  — 
n**  -  yj 

^  CloCvx 

cuCxx  = 


(15) 


w 


and  similarly  for  the  6  remaining  coefficients.  Here  C  is  the  radial 
bearing  clearance  and  W  is  the  load  on  the  bearing.  Introducing  Eqs .  (7) 
to  (10),  (14)  and  (15)  into  Eqs.  (12)  and  (13)  yields: 
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(16)  ["kjfx itoCxxJ  "h  C 1  Ijt  [Kjt^  ^  t loCky  ]  Ifo  -  H 

(17)  ^Kijy  +  UoCijX  ]  Xj  +  ~  Jd  +  ttoCjtj] 


where: 

azf 


The  bearing  housing  is  supported  flexibly  by  a  spring  Kp  and  a  dashpot  Cp. 
Denoting  the  mass  of  the  bearing  housing  hn  the  equations  of  motion  are: 

(19)  HiXj  =  ”  kp  Xj  ”  Cp  X'j  +  Ky* (Xj- X3 )  +  Cxlf  fxj  ”^3)  +  kxij  Cxif 


(18) 


X 


=  rc 

W<* 


Support  Equations 


and  similarly  for  the  y-direction.  To  make  dimensionless  set: 


(20) 

(21) 

(22) 


C  Co* 

w 


C<*>nC  P 

V 

-  Sll  J23L 
W*  f'M 


(23) 

(24) 


trt 


Solve  Eq. 
(25) 


(19)  together  with  Eq.(16)  and  (17)  to  get: 


*3  = 
^3  ~ 


jf(x*  +  0 


A-i  (g/+;(&w 

-  t) _ 


(26) 


System  Equations 


Combine  Eqs.  (16),  (17),  (25)  and  (26)  to  get  the  equations  in 
their  final  form: 

£  Kxx  ^  t  y  ( l+dxy )]  ^2  ^  itoCxij  tyj  = 

(27) 

£  Ito  Cijx  dijx  ]  Vj  "^  ["  Kt^  l  OjCtftj  ~~b(( I  1^2  ~  ~  ^  f  ^  l 


with  the  following  abbreviations: 


Kxx  “t  ItoCxx 


(28) 

§■ 

II 

k-  (&)*+  l  (&)cohc/ 

(29) 

k*4  ■+  i  Co  Cvif 

Oxij  ~ 

(30) 

kqx  ■+  t  Co  Cijx’ 

dyx 

Ki  +  '<»Cn  _  . 

(31) 

These  equations  may  be  solved  for  the  journal  amplitudes  x*  and  tj*  from  which 
the  rotor  amplitudes  and  the  bearing  amplitudes  can  be  determined  through 
Eqs.  (10) , (11) , (25)  and  (26).  The  solution  is  obtained  numerically  by  a 
computer. 

Singularities 


When  the  rotor  speed  CO  is  equal  to  one  of  the  rigid  support  rotor 
critical  speeds,  i.e.  ~  IjH  becomes  infinite.  If  6*X/  dxy  etc., 

are  finite  Eq.  (27)  yields: 


(32) 

X?  = 

(33) 

The  bearing  housing  amplitudes  are*found  through  Eq. 

I  +  <3^  •f  t  dxy 


(34) 


(35) 


*5  ~  '  "I"  (l  +  dxJCI+d^)-  4,^* 


^3 


*  if  I  - 


|-*-dxx  —  I  dlfx 


(l'fdx»)(l+dj(j)  —  dxydjx 


(19): 


X-  » 


The  rotor  mass  amplitudes  are  finally  given  through  Eqs .  (12)  and  (13): 


(36) 

(37) 


X,  =  -!  +f  [(Kw+Ud6o,)(*2-*3)+  (kMj  +  lCoC^Hifj -</,)] 

t],  =  t  •+  f  [(k^  +  iw<^J(xrx3)+flfy*  icuCW)(</i'-<ft)] 


‘  r*-y  (fenVM)Vo 


If  an  addition  to  y  =  00  >  <4Xv ,  <4xij  /  etc  are  also  infinite  the 
above  equations  reduce  to: 


(38) 

(39) 


x,  =  x2=x3  *-l 

Hi  '  Hi  "  Hi  ~  t* 


,  cond  =  0 


so  that  the  amplitude  is  180°  out  of  phase  with  the  unbalance  while 
the  center  of  gravity  of  the  rotor  does  not  move. 


Eqs.  (38)  and  (39)  are  actually  also  valid  even  if  as  l°ng  as 

k~  ,  tond~0  , 


Transmitted  Force 

The  force  transmitted  to  the  bearing  housing,  denoted  f  is  given  by: 

—  ^vvCXi~X3)  +  ( Xj  X3 )  "F  KXtj  Cxij  (^”^3) 

and  similarly  for  the  y-direction.  Making  use  of  Eqs. (12),  (13),  (10)  and 
(11)  this  equation  reduces  to: 

(4°)  ^  F8x  =  8e(l+x!) 

<«>  F’i  = 

The  force  transmitted  to  the  foundation,  denoted  F,  is  given  by: 
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which  through  Eqs .  (25)  and  (26)  becomes: 


(42) 

a  r  - 

Wcf 

(43) 

a  F  = 

W  6 

1  i 

1  o»d 

When  Xs  00  the  above  equations  become  indeterminant, 
use  the  equivalent  equations. 


In  that  case 


(“>  w  Fe<= 

(«) 

(44)  t>  =  [  4  +  /  v3 

(47>  wi  Fii  *  + 


Stability 

Since  the  support  flexibility  and  damping  influences  the  rotor-bearing 
stability  (oil  whip)  it  is  of  interest  to  analyze  this  effect.  To  do 
this  return  to  Eq.  (27)  and  investigate  the  eigenvalues  of  the  determinant 
of  the  homogenous  equations.  Introduce  the  symbols: 

(48)  jpw 

(49)  S*  & 


where  V  is  the  eigenfrequency .  Hence,  the  determinant  of  Eq.(27) 
becomes: 


(50) 


T C  +  toCxx  (1+6**)  ]  r  —  6xij  i 

[  ^.1  [  ktjtj  ■+  ~ J((l  + 


where: 


(51)  =  y- 

(52)  6*M  = 

(53)  6*h  = 


(sy)1 

Ms/ 


_ k»y  ~t~  t  Cxx _ 

h-  is%oo*d 

Kftf  ■^tyCoCvH _ 

4?-  ^f(Sj)l+  iSf(0*cl 


Look  first  at  the  case  where  the  supports  are  rigid,  i.e.  <4i)r  »  r  -- -•»(), 


Equate  the  real  and  imaginary  parts  of  Eq.  (50)  to  zero  to  get: 

_  Ko  Co  &L  4)C<y  Ky^  tp  CoCvij 

^  Cu  Cxx  -h£oC}) 


(54) 


(55) 


'i  CajCxxCoCui,^  OjC^u^/Cif^ 


Substituting  Eq.  (54)  into  Eq.(55)  allows  for  computing  y  ,  i.e.  the 
ratio  between  eigenfrequency  and  running  speed.  Substituting  into  Eq.(51) 
gives  the  value  of  S  j  i.e.  the  rotor  speed  at  instability  as  a  fraction 
of  the  natural  frequency  WM  . 


Next  assume  the  supports  to  be  flexible  but  without  damping,  i.e.  d-0 . 
For  abbreviation  set: 


(56)  A-  k-  ^(s/ 

and  the  solution  becomes: 


(57) 


(58) 


KxxCk)Cijij  ’j’HtfijUjCxY  kxij  (\)Cyx  ~~  k^r 

(ajCxx  4j Cijy 

{K^ZJElMaZ  T^F)- 

CoCxx  (oCljX 


i.e.  the  eigenfrequency  ratio  y  is  unchanged.  To  determine  the  rotor 
speed  at  instability,  i.e.  S  /  expand  Eq.  (57): 

(59) 

which  is  a  quadratic  equation  in  (5^)  .  Hence,  S 
y  known  from  EH.  (58). 


can  be  obtained  with 
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It  is  seen  that  there  are  two  solutions  for  the  speed  ratio  5  , 

both  with  the  same  frequency.  One  solution  corresponds  to  the 
coincidence  of  the  rotor  critical  speed  (including  the  effect  of 
the  flexible  support)  with  the  eigenfrequency  ^  ,  the  other  solution 
stems  from  the  instability  of  the  pedestal  mass.  This  was  checked 
numerically  as  follows: 

For  rigid  support  (  00  )  we  may  use  Eq.  (54)  and  (55)  to  determine: 

U>Q  J  '  COi, 

where  U>o  is  the  threshold  speed  for  rigid  support.  Let  COco  denote 
the  corresponding  actual  critical  speed.  Then.,  by  hypothesis: 


CO, 


'c,o  =  V 


i.t 


4 >c.o 

4;* 


-  M 


Introduce  an  equivalent  bearing  stiffness  Kg 


ils 

M 

Hence , 

V  Co*  1 

-  Md  -ft  -  c(  k. 

or 

(60) 

C  Kn 

W 

=  f  V 

such  that: 


Now  introduce  the  support  stiffness  to  get  an  effective  stiffness: 


(61) 


W 


ck,  _  3? 


ck  p 

w 


0<b  -tCKf 

w  w 


and  the  corresponding  critical  speed  becomes: 


*1  -  ft  -  v1 


Therefore . 


or 


/,  w)1  n  CKe 

'S,W  -  -  ?  W 


(62) 

which  is  the  first  solution.  To  evaluate  the  second  solution,  equate 
the  support  mass  resonance  to  the  eigenfrequency: 

+ 

h)  '  V 
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or; 


Thus: 


(63) 


CSr  +  Of a 
w  w 

hi 

M 


which  is  the  second  solution.  It  should  be  noted  that  both  Eqs .  (62) 
and  (63)  are  approximate.  Their  validity  has  been  checked  numerically 
for  selected  cases  and  the  difference  was  found  to  be  less  than  1  per 
cent..  Finally  when  the  support  includes  damping  set: 


(64) 


A-  k  -  ^  (ty)1  B-  (Sfr)cohc( 


Eq.  (27)  becomes 


(65) 


For  convenience  introduce: 


=o 


(66)  E-  ivF  -  r '”5T 

0  '  A-hB 

Expand  the  determinant  of  Eq.  (65)  in  its  real  and  imaginary  parts: 

(67)  F)(cuC^'i  F") ~ PCwij 

(68)  (Kxy“k  F)  +C  K|j^~E)(coC*x'*  F)  ky^  U>C  IjX  "*  kljy  ^oCy  ij  ~0 

Solve  Eq.  (68): 

_  kyy  (c<jC^'VF)~F  kij^  (  C«y  F  ^  Kyi|  ^Cjlf  4>Cyij 

(toCx,  -t  F )  + 


(69)  h. 
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Set: 

(70)  a  =  K^-K^ 

(71)  1 s  •=  CoC,* 

(72)  G|  =  Kyij  UjCij) (  ~f~  Kl j  y>  60  Cy  ^ 

(73)  X  -  F+  i(eoCkk-tCoC^) 


Substitute  Eq.  (69)  to  (73)  into  Eq.  (67): 

(74)  jV+  (Uw“4ofjj)]  =0 


This  equation  has  always  two  roots,  a  positive  and  a  negative  one.  Of  these 
only  the  positive  one  can  be  used  as  shown  later  in  Eq.  (79). 

Thus  Eq.  (74)  and,  therefore,  Eq.  (69)  can  be  solved  for  F  and  E  as  functions 
of  the  frequency  ratio  y  .  They  may  be  considered  as  frequency  dependent 
eigenvalues  such  that.  E  represents  "eigen  stiffness"  and  F  "eigen  damping". 
When  F  is  positive  the  "eigen  damping"  becomes  negative  according  to  Eq.(66). 
Having  determined  E  and  F  it  remains  to  find  the  corresponding  if -value 
from  Eq.  (66).  However,  it  should  first  be  noted  that  the  two  actual  un¬ 
known  are  ^  and  5  ~  coh  .  Since\the  calculation  is  performed  for  a 
given  Sommerfeld  Number  S  the  speed  Co  is  known  and  the  variable  is 
therefore  i.e.  the  rotor  mass  M  .  For  this  reason  it  is  convenient 

to  redefine  the  pedestal  mass  and  damping  as: 

4  =  f  s‘  -  “w1'  =  zn!5‘  (p^w) 

P  ‘  su"d  =  ^  (7057) 

Hence,  Eq.  (66)  may  be  written: 

(77)  E-ifrF  =• 

Solve  Eq.  (77): 


(75) 

(76) 


(78) 
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(79)  F"*‘r 

Combine  the  two  equations: 

<80,  K  -  Mjf^l 

Thus ,  for  a  range  of  frequency  ratios  ^  with  corresponding  E  and  F-values 
if  can  be  calculated  from  Eq.(80).  Substituting  into  the  right  hand  side 
of  Eq.  (79)  and  subtracting  from  F  gives  a  difference,  designated  the 
error.  Only  when  the  error  vanishes  has  a  solution  been  obtained,  i.e. 
for  a  correct  frequency  ratio. 

When  the  denominator  of  Eq.  (80)  has  a  singularity  a  false  solution  occurs 
except  in  the  rare  case  when  the  numerator  is  simultaneously  zero.  For 
this  purpose  it  is  convenient  to  have  a  check  on  Eq.  (80)  by  solving 
Eq.  (78): 

(8i)  +  £ -  0 


or 


<82>  *"  2[E^-<yn  [(*-<>»'> V+JE(H‘> 


Hence,  there  is  no  root  when: 


(83)  b  ~  i  E *  (^5)  *  E  ■+  /f1- 


Within  the  frequency  ratio  range  from  0  to 
be  up  to  three  roots  for 


frigid  ^Eq>  there  may 


CONCLUSIONS 


1.  A  flexible  bearing  support  provides  an  ideal  method  for  attenuating 
the  force  transmitted  by  an  unbalanced  rotor.  It  introduces  a  minimum 
of  extraneous  resonances  (namely  one,  the  resonance  of  the  bearing 
housing)  and  in  theory  the  attenuation  can  be  as  high  as  desired  by  a 
suitable  selection  of  the  support  stiffness. 

2.  The  damping  in  the  bearing  support  should  be  small  in  order  not  to 
lose  the  force  attenuation. 

3.  The  rotor  amplitude  at  resonance  is  relatively  high  but  since  it  occurs 
at  low  rotor  speeds  it  may  not  be  too  important.  Otherwise  the  support 
may  be  locked  in  passing  through  the  resonant  speeds. 

4.  Although  the  speed  at  onset  of  oil  whip  is  lowered  by  the  flexible 
support,  there  will  always,  in  practice,  be  enough  damping  in  the 
support  to  restore  the  stability. 


RECOMMENDATIONS 


1.  The  purpose  of  a  flexible  bearing  support  is  to  achieve  a  desired 
attenuation  of  the  transmitted  force.  However,  the  use  of  such  a 
support  raises  questions  with  regards  to  the  response  of  the  rotor 
during  onboard-ship  operation  (ship  motion,  shock  etc).  An  investi¬ 
gation  should  be  undertaken  of  this  problem  with  special  emphasis 

on  the  actual  design  of  the  flexible  support.  Important  factors  in 
the  study  should  be  reliability,  low  maintenance  and  simplicity. 

2.  Although  the  unbalance  is  the  principal  source  of  noise  in  a  piece 
of  rotating  machinery,  vibrations  of  other  frequencies  are  also 
present  and  may  be  important  (subharmonic  rotor  vibrations,  electric 
field  forces,  etc.).  A  study  of  these  vibrations  should  be  performed 
both  analytically  and  experimentally. 


-27- 


REFERENCES 

1.  Lund,  J.W.  and  Sternlicht,  B. ,  "Bearing  Attenuation,"  Bureau  of 
Ships  Report,  Contract  No.  Nobs  78930,  April  28,  1961. 

2.  Lund,  J.W.  and  Sternlicht,  B.  ,  "Rotor-Bearing  Dynamics  with  Emphasis 
on  Attenuation,"  Trans.  ASME,  Vol.  84,  series  D.  No. 4,  Journal  of 
Basic  Engineering,  Dec.  1962. 

3.  Warner,  P.C.,  and  Thoman,  R.J.,  "The  Effect  of  the  150-Degree  Partial 
Bearing  on  Rotor-Unbalance  Vibration,"  ASME  Paper  No.  63-LubS-6, 
presented  in  Boston,  June,  1963. 

4.  Warner,  P.C.,  "Static  and  Dynamic  Properties  of  Partial  Journal  Bearings," 
Trans.  ASME,  Vol:  85,  series  D.  No.  2,  Journal  of  Basic  Engineering,  June  1963. 

5.  Lund,  J.W. "Spring  and  Damping  Coefficients  for  the  Tilting-Pad  Journal 
Bearing,"  Bureau  of  Ships  Report,  Contract  No.  Nobs-86914,  March  1964. 


-28- 


APPENDIX 

Computer  Programs 

The  analysis  has  been  programmed  for  the  IBM  1620  digital  computer. 

Two  programs  are  written:  PN  0125  where  the  bearing  pedestals  are 
assumed  rigid,  and  PN  0132  where  the  pedestal'  flexibility  is  included. 
Separate  descriptions  of  the  programs  are  provided  below. 

PN  0125:  Transmitted  Force  and  Response  of  a  Two-Mass  Rotor  in 
_ Rigid  Pedestals _ 

This  program  calculates  the  dimensionless  transmitted  force  and  the 
rotor  amplitude  for  a  symmetrical  two  mass  rotor  with  a  given  un¬ 
balance.  The  rotor  is  supported  in  two  bearings,  each  bearing  re¬ 
presented  by  4  spring  and  4  damping  coefficients.  The  pedestals 
are  assumed  rigid. 


The  transmitted  force,  the  journal  amplitude  and  the  rotor  mass  ampli¬ 
tude  are  functions  of  the  rotor  speed  60  expressed  in  dimensionless 
form  by  (C%on)  where  0O„  is  the  critical  speed  of  the  rigidly 
supported  rotor.  However,  for  the  transmitted  force  and  the  journal 
amplitude  a  simplification  is  possible  by  use  of  another  variable: 


(a) 

where 

(b) 


*=7 


3  c? 

w 

c 

<* 

r 


1, 

-m 


(rotor  flexibilrtij 

bearing  reaction,  lbs. 
radial  clearance,  inch 


in 


60,  = 


influence  coefficient  for  rotor,  — s  (see  Eq.(4)) 

distance  between  rotor  mass,  divided  by  rotor  span 
(see  Eq.  (5)) 

critical  speed  of  rotor  on  rigid  supports, 
rad/ sec . 


2 

M  -  rotor  mass  per  bearing,  lbs-sec  /in. 


Since  is  inconvenient  to  use  directly  an  equivalent  speed  ratio 

may  be  defined  by: 

(c)  X=  10 

By  this  choice  of  non-dimensional  parameters,  a  single  plot  of  journal 
deflection  and  transmitted  force  against  equivalent  speed  ratio  will 
represent  a  full  range  of  rotor  flexibilities.  Since  the  rotor  de¬ 
flection  requires  the  designation  of  a  particular  flexibility,  the 
actual  speed  ratio  is  computed  in  this  case.  However,  for  large  values 
of  rotor  flexibility,  a  large  range  of  equivalent  speed  ratio  corresponds 
to  only  a  small  change  in  actual  speed  ratio.  Therefore,  in  this  case 
it  becomes  advisable  to  enter  the  actual  speed  ratio  and  so  this  pro¬ 
vision  is  made. 

The  program  also  performs  a  stability  analysis,  providing,  for  a  given 
rotor  and  bearing  configuration,  the  ratio  of  eigenfrequency  to  rotor 
speed  at  the  onset  of  instability  and  the  corresponding  rotor  speed. 

The  output  is  provided  in  a  form  that  makes  it  adaptable  for  automatic 
plotting  on  an  X-Y  plotter.  Since  it  is  anticipated  that  the  output 
would  be  represented  on  logarithmic  scales,  the  logarithm  of  the  re¬ 
sponse  values  along  with  a  code  for  automatic  sorting  is  provided. 

Computer  Input 

Card  1-49  columns  Hollerith  -  descriptive  text 
Card  2  -  (514,  E12.4)  -  control  parameters 

Word  1  -  type  of  computation;  if  this  value  rs: 

1.  the  input  speed  ratios  are  intended  to  be  equivalent 
speed  ratios 

2.  the  input  speed  ratios  are  intended  to  be  actual 
speed  ratios 

Word  2  -  number  of  speed  ratios,  maximum  18 

Word  3  -  number  of  rotor  flexibility  parameters  ,  maximum  10 

Word  4  -  intermediate  output;  if  this  value  is  0  -  no  inter¬ 
mediate  output  is  provided;  4  0  -  intermediate  output 
(diagnostic)  is  provided. 


(SLl 


Card  3 

Card  4 

Card  5 
Card 


Word  5  -  additional  data;  if  this  value  is  0  -  there  is  no 
additional  data;  4  0  -  additional  input  data 
follows  the  case  being  computed. 

Word  6  -  the  value  of  eccentricity  ratio  corresponding  to 
the  input  bearing  parameters. 


(6E12.4)  - 
efficients 

(6E12.4)  - 

(5E6.0)  - 

(6E12.)  - 


dimensionless  bearing  spring  and  damping  co- 
in  the  order  ,  <£*■ ,  ^  ^ 


the  remaining  bearing  parameters 


L y&L 

w 


CuiCfif 

W 


the  values  of  the  rotor  flexibility  parameter  Eq. 
the  speed  ratios 


Computer  output 

The  initial  output  is  in  the  order; 

-  program  heading 

-  Hollerith  text  provided  as  input 
input  control  parameter  values 

-  input  bearing  parameters 


Then  for  the  type  1  computation  the  first  output  values  with  appro 
priate  headings  are,  by  row 


first  row 


the  initial  equivalent  speed  ratio 

.  >vT 

the  maximum  dimensionless  journal  deflection 

the  value  of  the  dimensionless  parameter  3C  ,  eq.(&) 

a  column  of  signs  (+  1)  corresponding  to  the 
sign  of  the  log  of  the  speed  ratio  and  pro¬ 
vided  for  the  automatic  plotting  facility. 

the  log  of  the  speed  ratio 

the  sign  of  the  log  of  the  response  value 

the  log  of  the  response  value 

a  code  value  for  automatic  sorting 


second  row  -  the  same  as  the  first  row,  but  with  the  speed 
ratio  the  actual  speed  ratio,  the  deflection 
the  maximum  rotor  mass  deflection,  and  the  third 
value  the  value  of  rotor  flexibility  used  in  the 
computation. 

third  row  -  the  third  row  is  the  same  as  the  second  row  but 
for  a  second  value  of  rotor  flexibility,  if  more 
than  one  is  provided.  Then  the  additional  rows 
correspond  to  any  additional  flexibility  values. 


last  row  -  the  last  row  is  similar  to  the  above,  with  the  speed 
ratio  the  equivalent  speed  ratio  and  the  response 
value  the  dimensionless  transmitted  force  w 
(maximum  value). 

From  Equation  (a)  it  can  be  seen  that 

/W  _ 

l+*5> 


and  so  there  is  no  real  solution  to  for  values  of  iff  between  0 

and  -1.  For  those  cases  the  program  prints  the  message  "no  solution  for 
kappa  •  rho  between  0  and  -l  for  rho  =  .  Then  the  program  proceeds 

to  the  next  value  of  <j> . 


The  last  output  is  the  results  of  the  stability  analysis, 
of  three  values: 


It  consists 
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first  value  -  the  flexibility  parameter,^ 

second  value  -  the  speed  ratio  of  the  threshold  of  instability , 

third  value  -  the  square  of  the  ratio  of  eigenfrequency  to  running 
speed  at  the  threshold  of  instability,^} 


The  output  for  the  type  2  computation  is  much  like  that  of  the  type  1 
with  the  exception  that,  since  the  equivalent  speed  ratio  is  a  function 
of  both  the  actual  speed  ratio  and  the  flexibility  parameter,  each  time 
a  new  y -value  is  specified  for  a  given  actual  speed  ratio,  the  equivalent 
ratio  changes  and  requires  a  new  journal  deflection  and  transmitted  force 
computation.  Therefore,  the  order  of  output  is:  journal  deflection,  rotor 
deflection,  transmitted  force,  then  index  y  and  again,  journal  deflection, 
rotor  deflection,  transmitted  force,  etc. 


PN  0132:  Transmitted  Force  and  Response  of  a  Two-Mass  Rotor  in 
_ Flexible  Pedestal _ 

This  program  calculates  the  dimensionless  journal  deflection  and  trans¬ 
mitted  force  of  a  symmetrical  rotor  supported  by  two  similar  bearings  with 
flexible  pedestals  as  a  function  of  the  speed  ratio,  ^  ,  where  Coh  is 

the  flexural  critical  speed  of  the  rotor  simply-supported. 

The  program  also  performs  a  stability  analysis,  providing,  for  a  given 
rotor  and  bearing  configuration,  the  ratio  of  eigenfrequency  to  rotor 
speed  at  the  onset  of  instability  and  the  corresponding  rotor  speed. 

The  output  is  provided  in  a  form  that  makes  it  adaptable  for  automatic 
plotting  on  an  X-Y  plotter.  Since  it  is  anticipated  that  the  output 
would  be  represented  on  logarithmic  scales,  the  logarithm  of  the  re¬ 
sponse  along  with  a  code  for  automatic  sorting  is  provided. 


Computer  Input 


Card  1.-  49  columns  Hollerith  -  descriptive  text. 

Card  2.-  (614,  E12.4)  -  control  parameters 

Word  1  -  type  of  computation;  if  this  value  is: 


2  -  the  program  performs  the  rotor  response  computation  only 

3  -  the  program  performs  the  stability  analysis  only 

4  -  the  program  performs  both  analyses 
Word  2  -  the  number  of  speed  ratios,  maximum  18 


Word  3  -  the  number  of  non-dimensional  pedestal  to  rotor  mass 
ratios,  u,  where  n  =  J  maximum  10. 

Word  4  -  the  number  of  rotor  flexibility  parameters  y  ,  maximum  10 

Word  5  -  intermediate  output;  if  this  value  is  0  -  no  intermediate 
output  is  provided;  40-  intermediate  output  (diagnostic) 
is  provided. 


Word  6  -  additional  data,  if  this  value  is  0  -  there  is  no  add¬ 
itional  input  data;  4  0  -  additional  input  data  follows 
the  case  being  computed. 

Word  7  -  the  value  of  eccentricity  ratio  corresponding  to  the 
input  bearing  parameters. 


Card  3. 
Card  4. 


Card  5. 
Card 


-  (6E12.4)  *  Dimensionless  bearing  spring  and  damping  co¬ 
efficients  in  the  order  ^  yfe*  >  Sy  ' 

-  (6E12.4)  -  the  remaining  bearing  coefficients,  and  the  pedes¬ 
tal  parameters;  in  the  order 

where  the  parameters  are  dimensionless  as  defined  in  the  analysis. 

-  (5E6.0)  -  the  values  of  the  mass  ratio 

-  (5E6.0)  -  the  values  of  the  rotor  flexibility  parameter 


For  a  type  3  computation,  the  above  is  all  the  necessary  input;  other¬ 
wise  the  following  must  be  provided. 


Card  -  (6E12.4)  -  the  speed  ratios  (c3n) 


Computer  output 

The  initial  output  is  in  the  order: 
program  heading 

Hollerith  text  provided  as  input 
input  control  parameter  values 
input  bearing  and  pedestal  parameters 


Then  for  the  type  2  computation  the  next  output  with  appropriate  heading 

is  the  value  of  mass  ratio  and  rotor  flexibility,  followed  by  the  system 

response,  in  the  order 

-  first  row  -  the  initial  speed  ratio 

-  the  dimensionless  journal  deflection  IS  (max. value) 

d 


-  the  value  of  the  dimensionless  parameter 

.2 


where 


x  -A 


(-] 


*  I 


-  a  column  of  signs  (+  1)  corresponding  to  the  sign  of 

the  log  of  the  speed  ratio  and  provided  for  the  automatic 
plotting  facility 

-  the  log  of  the  speed  ratio 

-  the  sign  of  the  log  of  the  response  value 

-  the, log  of  the  response  value 

-  a  coding  value  for  automatic  sorting 


-  second  row 


-  the  same  as  the  first  row,  but  with  the  response 
value  now  the  dimensionless  rotor  deflection  ^ 
and  the  third  value  the  equivalent  speed  ratio. 

-  third  row  -  the  same  as  the  first  two  rows  but  for  the  di¬ 

mensionless  pedestal  deflection  £i»  (max  value) 

6 

-  fourth  row-  the  same  as  the  above  rows  but  for  the  dimensionless 

transmitted  force  I£E  (max.  values). 

W  6 

Then  the  next  speed  ratio  is  indexed  and  the  above  four  values  generated 
When  the  response  for  all  of  the  speed  ratios  are  computed,  a  new  rotor 
flexibility  value  is  indexed  and  the  above  sequence  repeated.  When  all 
of  the  flexibility  values  are  computed  a  new  mass  ratio  is  indexed  and 
the  above  process  repeated. 

For  the  type  3  computation  a  heading  indicating  the  stability  analysis 
follows  the  input.  Then  four  columns  of  output  as  follows: 

column  1  -  the  value  of  the  mass  ratio 

column  2  -  the  value  of  rotor  flexibility  factor 

column  3  -  the  ratio  of  running  speed  to  rotor  flexural  critical 

speed  at  the  threshold  of  instability 

column  4  -  the  square  of  the  ratio  of  eigenfrequency  to  running 
speed  at  the  onset  of  instability 

For  the  case  where  there  are  two  real  solutions  for  the  running  speed 
ratio,  they  are  both  given. 

For  the  type  4  computation  the  output  after  the  listing  of  the  input 
values  is  the  type  2  output  followed  by  the  type  3  output. 


Dimensionless  Support  Stiffness  =  1C" 
Dimensionless  Bearing  Mass  =  0.C5 
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SPEED  RATIO  w  Pig.  4  Dimensionless  Transmitted  F  ,r  :e  vs  Rot  ,r  Speed  Ratio 

Dimensionless  Support  Stiffness  =  3  ■  ic_3 
Dimensionless  Bearing  Mass  =  C.''5 


SPEED  RATIO  -g-  Fig.  5  Dimensionless  Transmitted  Force  vs  Rotor  Speed  Ratio 
wn  Dimensionless  Support  Stiffness  -  3  •  i:*3 

Dimensionless  Bearing  Mass  =  0.15 


Dimensionless  Support  Stiffness  =  10 
Dimensionless  Bearing  Mass  =  0-C5 


Fig.  7  Dimensionless  Transmitted  Force  vs  Rotor  Speed  Ratio 
Dimensionless  Support  Stiffness  =  10~2 
Dimensionless  Bearing  Mass  =  0.15 
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Dimensionless  Transmitted  Force  vs  Rotor  Speed  Ratio 
Dimensionless  Support  Stiffness  =  3  •  10'2 
Dimensionless  Bearing  Mass  =  0.G5 


litted  Force  vs  Rotor  Speed  Ratio 
•t  Stiffness  =  3  •  10~2 
ig  Mass  =  0. 15 
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Fig.  12  Dimensionless  Journal  Amplitude  vs  Rotor  Speed  Ratio 
Dimensionless  Support  Stiffness  =  3  -  10"3 
Dimensionless  Bearing  Mass  =  0.C5 
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pig.  23  Dimensionless  Rotor  Mass  Amplitude  vs 
Dimensionless  Support  Stiffness  =  10"' 
Dimensionless  Bearing  Mass  =  0. 15 
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NOHENCLATURE 

C 

Cxx  jCxij/  ClJX,  Clflf 

&>C(X 

Cp 

D 

cohd 

E 

F 

F 

Ft  /  F 

Kor ,  K^,Kqx, 

Kw 


'XX 


kP 

A 

L 

/ 

M 

N 

R 

S 

5 

t 


Radial  bearing  clearance,  inch 

Bearing  damping  coefficients,  lbs-sec/inch 
_  C  Co  Cxx 

=  t  dimensionless  bearing  damping 

coefficient.  Similar  for  &>Q(y,CuGjy,£oG|ij 

Pedestal  damping  coefficient,  lbs-sec/inch 

Bearing  diameter,  inch 

-  C(OnCf 


W 


t  dimensionless  pedestal  damping 
coefficient 


Real  part  of  eigenvalue  of  homogenous  equations 
Imaginary  part  of  eigenvalue  of  homogenous  equations 
Transmitted  force,  lbs. 

V  AhJ  (j- components  of  transmitted  force,  lbs. 

Bearing  spring  coefficients,  lbs/in. 

*  Cw“  dimensionless  bearing  spring 

coefficient.  Similar  for 

Pedestal  spring  coefficient,  lbs/in. 

-  ^dimensionless  pedestal  spring  coefficient 

Bearing  Length,  inch 

Rotor  span  between  bearing  centers,  inch 

2 

Half  the  rotor  mass,  lbs-sec  /in. 

2 

Mass  of  bearing  housing,  lbs-sec  /in. 

Rotational  speed  of  rotor,  RPS 

Bearing  radius ,  inch 

ixNDL  /R\2 

~  y/  \C  /  /  Sommerfeld  number 
£0 

s  ton  )  Ratio  of  speed  at  instability  and  rotor 

critical  speed 

* 

Time,  seconds 


-si- 


vy 

Xi/ij i 

Xi/Hi 

*3/  Hj 
d&a.  )  ^afc 


Bearing  reaction,  lbs. 

Amplitude  of  rotor  mass,  inch 

Amplitude  of  journal  center,  inch 

Amplitude  of  bearing  housing,  inch 

”  ^  )  Dimensionless  amplitude  of  rotor  mass 

Dimensionless  nirpiltade  of  journal  center 

P;  Dimensionless  amplitude  of  bearing  housing. 

Rotor  influence  coefficient,  first  index:  amplitude, 
second  index:  force,  in/lbs. 


* 

/» 

Jf 

6 

x 

h 

/* 
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*  (oUa+oUi)  for  static  unbalance,  ■  Co<aa— oCat)  for 
dynamic  unbalance,  in/lbs. 

_  CmCp 


W 


Dimensionless  pedestal  damping  coefficient 


=  ^  Ratio  of  instability  frequency  and  rotor  speed 
Unbalance  eccentricity  of  rotor  masses,  inch 


parameter 


Lubricant  viscosity,  lbs-sec/in 

rr\ 

=  fFi  ) 

Journal  eccentricity  ratio  in  bearing 


Dimensionless  bearing  housing  mass 


V 


Instability  frequency,  rad/sec 


! 

? 
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*  1  for  static  unbalance,  «*  ratio  of  distance  between 
rotor  masses  and  total  rotor  length  l  for 
dynamic  unbalance. 


Wot 


w  > 


rotor  flexibility  parameter 
dimensionless  bearing  housing  mass 


dxx i^xi] /dijg'dyy  Parameters  defined  by  Eqs.  (28)  to  (31) 


(0 


Angular  speed  of  rotor,  rad/sec. 


critical  speed  of  rotor  on  rigid  supports, 
rad/sec. 
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SUMMARY 


The  dynamic  response  of  an  unbalanced  elastic  rotor  which  operates  in  damped 
fluid-film  bearings  has  been  investigated.  The  influence  of  speed,  bearing 
operating  eccentricity,  and  of  relative  stiffness  between  the  rotor  and  its 
bearings  have  been  determined.  Results  are  expressed  in  terms  of  rotor  maximum 
whirl  amplitude  at  several  locations  on  the  rotor;  and  in  terms  of  bearing 
transmitted  force.  Particular  attention  has  been  given  to  the  influence  of 
system  parameters  on  critical  speeds,  and  on  the  attenuation  of  bearing  trans¬ 
mitted  force. 

The  results  are  presented  as  charts  which  facilitate  the  design  of  high-speed 
rotors  by  allowing  performance  characteristics  up  to  and  including  the  fourth 
system  critical  to  be  determined  directly.  The  analysis  is  exact,  and  the  influ¬ 
ence  of  higher  modes  is  included  under  all  operating  conditions.  Both  static 
and  dynamic  unbalance  conditions  have  been  considered. 
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INTRODUCTION 

General 

A  heavy  elastic  rotor  has  an  infinite  number  of  critical  speeds,  but,  in  practice, 
the  operating  speed  range  rarely  includes  more  than  three  or  four  of  these  criticals. 
These  critical  speeds  and  the  vibratory  form  assumed  by  the  rotor  are  both  directly 
influenced  by  the  stiffness  and  damping  properties  of  the  supporting  structure, 
as  well  as  by  the  distribution  of  mass,  elasticity,  unbalance  and  damping  within 
the  rotor  itself. 

Little  information  is  presently  available  on  the  dynamic  response  of  rotor 
systems  to  unbalance  in  the  higher  modes.  The  major  studies  which  have  been 
made  to  date  have  been  concerned  with  optimising  the  rotor-bearing  system  to 
achieve  maximum  attenuation  of  transmitted  force  and  rotor  amplitude  in  the 
lower  modes.  Lund  and  Sternlicht  (Ref.  1)  used  a  simple  rotor  in  fluid-film  bearings 
to  investigate  force  attenuation  in  the  fundamental  mode  for  a  variety  of  bearing 
types.  Warner  and  Thoman  (Ref.  2)  extended  this  work  using  a  two-mass  rotor  in 
partial-arc  bearings.  This  rotor  included  the  influence  of  the  second  mode  in 
the  fundamental  rotor  motions,  and,  also,  provided  data  on  operation  at  speeds 
including  the  second  critical. 

In  both  cases,  the  simplicity  of  the  rotor  precludes  any  direct  extension  of  this 
work  to  higher  modes.  However,  such  information  is  desirable  in  order  to  opti¬ 
mise  the  steady-state  performance  of  high-speed  rotors  which  operate  beyond  the 
second  bending  critical,  and  to  optimise  the  system  run-up  and  run-down  charac¬ 
teristics.  Also,  the  influence  of  the  higher  modes  on  rotor  response  in  the  lower 
modes  and  throughout  the  speed  range  is  not  indicated  in  presently  available  data. 
Such  information  is  required  to  enable  the  most  efficient  attenuation  to  be  deter¬ 
mined.  This  information  is  obtained  in  the  present  analysis.  Finally,  the  para¬ 
metric  conditions  under  which  stable  rotor  motions  may  occur  during  operation  in 
the  higher  modes  have  not  been  defined  up  to  the  present.  A  simple  extension  of 
the  results  given  herein  would  allow  this  to  be  done. 

Scope  of  Present  Investigation 

In  the  present  analysis,  the  motions  of  an  unbalanced  flexible  rotor  which  is 
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supported  in  fluid-film  bearings  are  considered.  The  rotor  is  assumed  to  be  of 
prismatic  shape,  and  to  have  its  mass  and  elasticity  distributed  uniformly  along 
its  length.  This  distribution  makes  it  possible  for  the  solutions  to  the  result¬ 
ing  equations  of  motion  of  the  system  to  include  the  influence  of  all  modes 
directly.  Rotor  internal  damping  is  assumed  to  be  negligible  compared  with  the 
damping  in  the  fluid-film  bearings.  A  simple  uniform  shaft  was  chosen  because 
it  permits  a  direct  analytical  solution  to  be  obtained  for  all  required  dynamic 
properties  of  the  system  including  the  higher  modes.  This  rotor-bearing  system 
is  not  intended  to  simulate  any  practical  case,  but  it  does  indicate  certain 
dynamic  characteristics  which  are  common  to  all  cases. 

The  motions  of  the  rotor  are  considered  to  arise  from  the  action  of  an  unbalance 

W  *  e  which  is  due  to  the  weight  of  the  rotor,  W,  acting  at  an  eccentricity  e, 

from  the  rotor  geometric  centerline.  The  unbalance  is  located  at  a  specified 
point  along  the  length  of  the  rotor,  and  its  distance  from  the  left-hand  bearing 

is  included  as  a  variable  in  the  analysis.  This  is  a  "static"  unbalance.  Any 

desired  condition  of  dynamic  unbalance  may  also  be  investigated  by  suitably 
superposing  two  sets  of  static  unbalance  results  to  represent  the  dynamic  un¬ 
balance  couple.  This  superposition  has  been  included  in  the  analysis,  and 
dynamic  unbalance  results  have  been  obtained.  The  point  at  which  the  desired 
amplitude  or  force  response  occurs  is  also  included  as  a  variable. 

The  rotor  is  supported  in  a  hydrodynamic  cylindrical  fluid-film  journal  bearing 
at  either  end.  The  dynamic  properties  of  these  bearings  are  given  in  Table  I. 
This  bearing  type  is  commonly  used  in  practice.  It  is  considered  in  the  ana¬ 
lysis  by  Lund  and  Sternlicht  (Ref.  1),  and  its  characteristics  are  similar  to 
those  of  the  partial-arc  bearing.  Thus,  this  choice  of  bearing  allows  a  com¬ 
parison  to  be  made  between  the  present  and  previous  work. 

The  analysis  presented  here  has  been  programmed  and  the  following  results  have 
been  obtained: 

1.  a.  Rotor  amplitude  at  specified  stations  for  a  speed  range  which  in¬ 

cludes  the  fourth  rigid  bearing  critical, 
b.  Bearing  transmitted  force  for  the  same  speed  range. 

2.  Rotor  mode  shape  in  terms  of  an  optional  number  of  rotor  stations  for 
selected  speeds. 
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DISCUSSION 

Non-Dimensional  Parameters 

All  parameters  in  the  results,  Figures  6  through  24,  have  been  made  non-dimen¬ 
sional  for  generality.  The  parameters  and  dimensionless  ratios  used  are  de¬ 
fined  as  follows; 


1.  Notation.  Listed  at  the  end  of  this  report 

2.  Distance  ratio,  L^/L  =  aC  .  Defines  the  axial  location  of  the  unbalance 
W  •  e  from  the  L.H.  bearing,  with  respect  to  shaft  length.  Figure  1. 

3.  Position  ratio,  Z^/L^.  Defines  the  axial  position  of  any  displacement 
measurement  within  Region  1,  measured  from  the  L.H.  bearing.  Figure  1. 

4.  Position  ratio,  Defines  the  axial  position  of  any  displacement 

measurement  within  Region  2,  measured  from  the  unbalance  position. 

Figure  1. 

5.  Speed  ratio,  .  Ratio  of  rotor  speed  co  to  the  fundamental  bending 

critical  speed  cdc  of  a  uniform  rotor  in  simple  rigid  end  supports.  The 
speed  range  covered  by  the  analysis  extends  to  oj/o>c  *  24.0.  This  includes 
the  first  four  rigid  bearing  criticals. 

For  the  present  rotor  system  co  = 


where  \L  is  the  system  characteristic  frequency  number. 

3t2  r  Ell  ^ 

For  a  rigid  bearing  uniform  rotor  co  =  — r  — 7 

c  L2  IpAJ 

Speed  ratio  <c~)  =  (~)^  At  the  rigid  bearing  critical  co/co  =  1.0. 


Characteristic  number  \L  =  n  (—) 

CJD 

c 


A.L.  =  £  it  (^-)‘ 

1  <0 


\l2  =  (i-JC)it  <-*-)' 


and 


-4- 


6.  Flexibility  ratio,  — .  Ratio  of  the  bearing  clearance  C  to  the  central 

0 

deflection  of  a  uniformly  loaded  simply- supported  shaft,  5. 

3  C 

7.  Stiffness  parameter,  v  *  EIX.  — .  Occurs  non-dimensionally  in  the  equa- 

"  Q 

tions  of  motion.  The  relationship  between  V  and  —  is  as  follows: 


Central  deflection  of  simply- supported  uniform  shaft  5  = 


Static  load  per  bearing 
Characteristic  number 
3 


W, 


5  VAi 
384  El 

WAL 

2 


XL  -  «  (B-)4 
'u> 


Therefore  v  =  EI\  • 


3  3/2 

C  .  (52L.)(C)(cl;^ 

W  '192  'V  to  • 
c 


Stiffness  parameter  is  therefore  a  dynamic  deflection  ratio,  as  shaft 
stiffness  depends  on  rotor  speed.  The  range  of  flexibility  ratio  re¬ 
quired  for  the  analysis  was  established  as  follows: 

Bearing  Fundamental  Rotor  Shaft  Lateral  C 

Case  Clearance  C.  in  Critical  Speed,  rpm  Deflection. in.  8 

Maximum  0.0050  14,000  0.002  30 

Minimum  0.0005  4,000  0.015  0.3 


The  fundamental  rotor  critical  speeds  were  chosen  high  in  order  to  apply 

to  a  simply- supported  shaft  in  rigid  bearings.  Flexible  bearing  system 

c 

criticals  will  be  considerably  lower.  The  ~  values  therefore  will  apply 

0 

to  the  range  for  rotor-bearing  system  fundamental  critical  speeds.  The 
£ 

above  maximum  —  value  therefore  applies  to  a  fairly  rigid  shaft  in  flex- 

8 

ible  bearings.  The  lowest  critical  speed  will  tend  towards  a  rigid-body 

critical  and  will  be  determined  mainly  by  bearing  flexibility  and  rotor 

C 

mass.  The  minimum  value  —  =  0.3  corresponds  to  .a  flexible  shaft  in 

5 

rigid  bearings  The  lowest  critical  speed  will  tend  to  occur  at  oj/co^  =  1.0, 
determined  by  shaft  flexibility.  These  facts  are  of  use  in  analyzing  the 
results  given  in  Figures  6  through  24 


8.  Displacement  amplitude  ratio,  — .  Ratio  of  the  maximum  rotor  displacement 
x  at  a  specified  station  to  the  unbalance  eccentricity  e. 

•CF 

9.  Dimensionless  transmitted  force,  77-.  Ratio  of  the  maximum  transmitted  force 

We 

F  to  the  unbalance  W  •  e,  normalized  by  the  inclusion  of  the  bearing  clear¬ 
ance  C. 
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Bearlng  Properties 

The  bearing  used  in  the  present  analysis  is  a  cylindrical  journal  bearing. 
Dynamic  stiffness  and  damping  characteristics  are  given  below  in  Table  I, 
derived  from  Reference  1  in  Appendix  A  of  this  report. 


TABLE  I 

Spring  and  Damping  Properties  for  Cylindrical  Journal  Bearing,  L/D  =  1 


0.2 

0.5 

0.7 

K 

XX 

1.283 

2.060 

3.59 

K 

xy 

5.492 

3.230 

3.38 

K 

yx 

-4.610 

-1.070 

0.02 

K 

yy 

2. 220 

2.040 

1.99 

c 

XX 

10.72 

6.02 

6.23 

C 

xy 

1.950 

2.00 

1.95 

C 

yx 

2  290 

2.17 

2.13 

G 

yy 

9.770 

3.40 

2.00 

s 

0.665 

0.189 

0.081 

76. 8 

57.5 

43.3 

The  cylindrical  journal  bearing  is  similar  to  the  partial-arc  bearing  in  opera¬ 
tion  because  of  the  presence  of  cavitation  in  the  film.  Thus,  the  results  ob¬ 
tained  for  the  system  characteristics  are  representative  for  a  wide  range  of 
common  applications.  For  systems  which  employ  bearing  types  with  important 
differences,  such  as  a  tilting-pad  bearing,  the  qualitative  aspects  of  the  results 
still  apply  as  a  guide  even  though  they  are  not  then  numerically  correct. 

Maximum  Response  Values 

Response  amplitude  and  transmitted  force  results  are  both  given  for  the  maximum 
values  for  these  parameters.  The  method  of  calculation  for  both  maxima  types 
is  given  in  the  theoretical  analysis,  Section  3,  and  a  diagrammatic  representa¬ 
tion  of  the  component  relationships  which  go  to  make  up  the  force  and  displace¬ 
ment  maximum  values  is  given  in  Figure  25 
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Both  static  and  dynamic  unbalance  characteristics  are  included  in  the  above  re¬ 
sults.  Details  for  Program  1  are  given  in  Figures  6  through  24,  and  for  Program 
2  in  Figures  26  and  27.  All  major  system  variables  are  expressed  in  terms 
of  dimensionless  parameters 

The  results  of  the  present  investigation  have  been  compared  with  those  obtained 
using  a  proven  discrete -mass  rotor-bearing  program.  Correlation  was  good 
in  all  cases.  The  results  are  shown  in  Figure  28. 


-7- 


THEORETICAL  ANALYSIS  OF  ROTOR  MOTIONS 
General 

The  rotor-bearing  system  shown  in  Fig.  1  consists  of  a  heavy  elastic  shaft  of 
uniform  circular  cross-section  along  its  length,  supported  at  its  ends  in 
fluid-film  bearings.  The  bearings  have  stiffness  and  damping  properties  in  both 
x-  and  y-  directions,  and  also  cross-coupling  stiffness  and  damping  between  these 
directions.  Dimensionless  values  of  the  bearing  stiffness  and  damping  properties 
with  eccentricity  are  given  in  Table  1  for  a  cylindrical  journal  bearing. 

The  shaft  considered  has  an  unbalance  W-e  situated  at  distance  L1  from  Z,  *  0. 

2  i  1 

During  rotation,  this  gives  rise  to  an  unbalance  force  Meco  which  rotates  in 
synchronism  with  the  shaft,  causing  it  to  whirl  about  its  stationary  equilibrium 
position.  Shaft  motions  are  restrained  by  its  own  inertia  and  elasticity,  dis¬ 
tributed  uniformly  along  its  length;  and  by  the  bearings  at  either  end.  The 
bearing  fluid  film  forces  consist  of  a  linear  spring  force  which  opposed  dis¬ 
placements,  and  a  viscous  damping  force  due  to  velocity.  Any  externally  impressed 
journal  motion  in  a  given  direction  gives  rise  to  fluid-film  forces  which  oppose 
the  motion,  both  in  the  direction  of  the  displacement  and  at  right  angles  to  it. 

The  coordinate  bearing  forces  arising  from  journal  motion  are  written  as 

y-d.  rection:  F(Y)  =  Kyy  Y  Cyy  Y  +  K*y  X  ♦  C*yX 

x-direction:  F(*)  *  Kxx  X  ■+  *■  Ky*  Y  +  Cy*Y 

The  unbalance  force  may  also  be  resolved  into  the  x-  and  y-  directions.  This 
allows  the  usual  equations  for  plane  motions  to  be  written  for  rotor  displacements. 
The  solutions  to  these  equations  may  be  combined  to  yield  the  maximum  displacements 
and  forces  acting  on  the  rotor. 

The  purpose  of  this  investigation  is  to  examine  the  nature  of  shaft  displacements 
and  transmitted  bearing  forces  which  result  from  shaft  unbalance  in  the  system 
described,  over  a  wide  range  of  speed. 
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Baslc  Equations  and  Solutions 


Considering  motion  in  the  y-direction  of  Figure  (1),  the  equilibrium  of  an  elemental 
length  of  shaft  dz  subject  only  to  internal  forces  is  governed  by  the  well-known 
equation 


\\ 

eiH4  + 

•urA 

1 

0 

i 

>/-p 

/o'''® 

(i) 

The  solution  to  Equation  (1)  is 

T 

,  N 

-  ^(*>e 

(2) 

where  y(z)  is  a  function  of  z  and  independent  of  t.  For  the  section  0(2,  <  L, 
substitution  gives 


-  -i  ixwt  ( 3 ) 

Yt  S  I  AjCosXz*  ♦  +  C,co*hX-*\  *  0,  SmhXz,  J  e 

.4  vA  <ox 

*  *  9-ei 

and  Dj  are  constants  of  integration  to  be  determined  from  the  end 
the  section.  For  the  section  0  i  3^  <  Lj.  the  solution  is 

*  [Ai&stai*  6*5* Azx  +  Cj.cosh*zx  +  ^swhAZi]  £*  (4) 

C 2  and  are  constants  of  integration. 


where 

and  A^, 
conditions  of 

where  A2>  B2, 


Similarly,  for  motion  in  the  x-direction,  the  governing  equation  is 


For 


O  ^  Z,  *  L 


eisX  ♦  -  o 

iz4  ^  it 

the  solution  is 


(5) 


X1,  *  [  E,  «0  Az,  <-  Sin  AZ,  «■  <*t  C»Sll  AZ,  4  H,  SwlA  A*,]  e 


hot 


(6) 


where  E^,  F^  and  are  constants  of  integration  and  \  is  as  defined  above. 
For  O  4  Zx  *  1-2.  the  solution  is 

Xz  »  J^ElcosAzl+  <sxcoshAzL  f  Mj.  SnAAZjj  &**** 

where  E2>  F2>  G2>  and  H2  are  constants  of  integration. 
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Evaluatlon  of  Constants 


The  16  constants  of  integration  may  be  evaluated  by  introducing  the  solutions  con¬ 
taining  them  into  the  boundary  conditions  of  the  system.  Adopting  the  conventions 
of  Figure  2  allows  the  boundary  conditions  to  be  expressed  as  follows: 


i  y-z  f* 


fane 


KC'OCo')  *  o  (8) 

V(V,)(o)  -  o  (9) 

where  F(XXo)  =  Kyy  X  (<0  *  Y,(°)  +  KxyX,(o)  +  X.C©)  (10) 

21  *  °  i  I**"*  M(X,Xc>)  -  O  (U) 

vCx.Xo)  *  o  (12) 

where  F  =  Kot^tC0')^  C*  X,( o)  +  K^XCp)  +  Cy*  C°)  (13) 

1-1  J  y-Z  plane.  mCyjco  -  o  (H) 

mooo  +  fcyoc«-o  *  o  (is) 

where  F(XX*a)  *  Kyy  Yx(Ll')  +  Cyy  Y2  (Lv)  +  ^  Xx0*)  +  Qy  Xj(Lx)(16) 

Zi*  Ll  >  aC-Z  plane.  M  (jQCU)  *  O  (17) 

vcxjcu.)  +  FCOO-0  *  °  (18) 

F(XXW  -  kM xt00  .  +  KwXc*-*)  *  cwt(0<19) 


where 


Z,  *  L,  :  ij-z  plane. 
2^0: 


r<»-) 


X  (o) 


(20) 

(21) 

(22) 


im  - 

mCy.Xli)  =  m(xX°> 

-  VCY.XO  ♦  Hewv[-ieiwfc]  4  Vfc/o).  o  (23) 

where  Mcu?"  [-16  J  is  the  y-  component  of  the  unbalance  force  acting  at  . 

x,M  -  XtCo) 

=  4.oc‘-x°')  ta) 


• 

t  e 

2(*L,  :  a-z  pane. 

Z*  *  0  j 
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MCXiX  L0 

«  M(X*Xo) 

$0 

-V(X,XO  +  Mew  •  emt  4 •  VfoX6) 

-  0 

<&) 

where  MetO*^**  is  the  x-  component  of  the  unbalance 

force  acting  at  L^. 

Substituting  eqp.  (3)-(7)  int  eqs.  (8)-(27),  and  utilizing  the  moment  and  shear 
relations 


and 


m(y)  - 

dV 

'  er 

M(X) 

VCY)  * 

dS 

E1  dz3 

vOO 

-  El 


dor 

di'Z*- 


-  El 


cLz3 


allows  the  16  basic  equations  of  motion  to  be  obtained. 

Expressed  in  dimensionless  form  in  terms  of  the  integration  constants,  these  equa¬ 
tions  are  as  follows: 


A, 

e 

9 

e 

*  O 

«.($  ♦§') 

— 

i>  (|  -  | )  *  *»  ( t  ♦  t ) 

=  o 

Cl  -  !') 

»  o 

4 

^4  C|  ♦ 1)  -  ■»  Cl  -  1) 

=  o 

^  CoS  >l-i 

+ 

^  Sin  —  £?■  cosh  ^  Swh  At*  *  O 

e  «  e 

-1  4 

e  5 

1  **  *  t  *>  ♦  |  *« 

+  <*i 

[1 

cwAij.  4  5  Aii  4  cosh  k.x  4 

Sinh  Ati  j  =  O 

1* 

+ 

§  sm  ALi  -  <!?■  cohAl^  -  ^ 

e  e  e 

SinA  Atj.  *  O 

[ 

'  Ai 
•  « 

Co  AL*  +  ^  «n  A^  +  ^  c»#h.  At*  4 

^  Sinh  1 

*■  %  *»  - 

'  *.o  +  ^  4  £ 

°<n.  =  O 

6  »AU, 

♦ 

§  Sin  At,  4  ^  c*oh,  At,  4  ^ 

c  C  g 

SmA  >t, 

A> 

e 


Cl 


O 
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Ai 

e 


?!  a* 

c 


4.  Si nk  Ai-i 

e 

-  b 

e 


+  B  c&sM  M-i 
e 


—  OA  Xi 

Si*  AL| 

|  asAf-, 
e 

-I  «.*, 

-  b  & s 
e 

£•  s/o  A*-i 
6 


-  i>  SIM  Al, 

e 

+  -1  CoSh  A*-| 

4.  2f  smh  A^i 

£» 

+ 

e. 

C(  «  O 

e 

-  §  cos  Al, 

c 

♦  |  Si*AAL, 

♦  1 

4.  B  cosh  At-i 
e 

9»  =  tAZ- 

c 

4-  ^  Sim  A*-/ 

w 

+  &  Cosh  Mi 

e 

4-  (ff  S#»A  A*-! 

45 

-  £ 
e 

<5x  *  O 

e 

4-  5  COS  V/ 
£ 

+  Sinh  Mi 

e 

-  £  - 
e 

4-  COS^J  M-i 

e 

*2.  =  O 

e 

-  a  sw  M-i 
e 

+■  -'  cuih  Mi 

e 

+  /b  swA  AZ-/ 
e 

*  ei 

<*2.  sz  O 

e 

-  5  CoS  M,  +  t'  *">* 
e  * 

+  b  ctiih  Mi 

e 

+  I 

e 

_  ^  »  —  az. 

e 

(28) 

*■(«!)*, 

iasst  W 

(EIf«  J 

A*-i  = 

Jf  *  <£  >4 

y 

r 

II 

Mi*  (ft)4 

•j)  = 

(gxi)(V* 

EZ>?C 

W 

in  which 
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<*\ 

,  r  *&Jr  +  l. 
1  w 

CcoCyyl 

W  J 

-  [<2Lf  +  i. 

L  W 

CwCcy] 
w  J 

*3 

=  1*  k'*c  i 

L  w 

CtoCyi  1 
W  -1 

°<4- 

*  1"  *£L?  +  i 

.  Cw  C**  1 
W  -1 

*  £  -  V  sm  >/.a 

+  o^|  COS  )li-j  J 

°<6 

=  [  V  c»S  >LX 

+  <*,  sir,  >LX] 

<*7 

«  [-?)  sw#>  )LX 

■*•  o^|  66S#,  ^Z.  x  J 

*8 

-  [-V  COS#,  ^ 

+  o^|  S,n#»  AZ-x^J 

*  [-V  "Sin  )LX 

+  o<4 

COS  ^Z-x  ] 

<*to 

*  [  V  cos  >LX 

+  °<4 

Sin  >4 

°<ll 

»  [-7>  &nh  kx 

+  *4 

cz**'  AZ-x  ] 

*11 

=  cosh  k: 

l  +  °<4 

Sin#i  Mx  J 
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It  Is  convenient  to  reduce  the  preceding  equations  to  a  set  of  6  simultaneous 


equations  In  ^  E^  to  facilitate  the  computer  solution  which  Is  to 

follow.  After  reducing  and  simplifying  the  equations  become: 

-  |-v 

4- 

B.V 

e 

4-  £'.2*!  *  O 

%-2«)  * 

— 

5.v 

e 

4.  H«.v  »  O 

^  [  CoS  Ad  -  toh  Ad  J  4  |  smAl 

— 

B  Smh  AL 
e 

4  L  (%)[*"»  Ad*  +  »«AAdJ 
*  O 

h  [cosAl  -  Cosh  Ad]  +  §  siflAi- 

— 

M*  SnA  Ad 
€ 

—  (4-j  )  [sn  Ada  4  AlxJ 

=  O 

[o<y  OtfAd,  -  +  oC7  C0S#»  \L, 

4.  [oC7  swh  Ad,  4  <x6CoshH,]  + 

+  o^WiAd/J  +  |  [<xsswAd,  *  o^cosAz.,J 
&  o(x  [d»S  Ad  cash  Ad.  j  4-  ^o^SihAd  4  Mp  SmAA 

*  (£)  [o^(ah 


&  [swAC  ♦  Ad  J  +  §  <*a  Swi  At  4  B  ^  sJnh  A/- 
+  j|  [o^cosAd*  +  «X/i  £®s/»  Ad./  +  od/Jt  S*AAd/J  4  £  [ 

+  *  fasinhk,  4  Of^coihk,]  *  (£  )[«-*»)-  £*b(sh *2 “**>*■')]  (29) 

The  following  expressions  exist  for  the  remaining  constants  in  terms  of  the  above 
unknowns . 

A,  C, 

e  "  c 

i»  =  £. 

e  e- 
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tg.  3  CM  XL,  4-  |  SinXL, 

^  s  -  6  sin  XL,  4  S  OS^i  +  i  (  ^  ) 

—  =  Ai  cosh  XL,  4  B  swhAL, 

e.  £  e 

^  a  £'  Sinh  XL,  4  B  cO<h  XL,  -  -t  ( 

^  ^  coS  )L|  +  5  Sm  XL, 

e-  e  2 

^  a  -  §  an>,  4  5  COS  XL,  -  (^) 

*  e  e 

a  &  cosh  XL,  4  *  Smh  XL| 

C  €  6 

y?  -  §1  smh  >L,  4*  Hi  cosh  XL,  4-  } 

Ci  6i  £ 


(30) 


Calculation  of  Rotor  Displacements 

In  the  region  O  S  ^  S  L(  ,  the  respective  coordinate  displacements  are 
given  by  Eqs.  (3)  and  (6)  In  order  to  determine  the  maximum  and  minimum  values 
of  whirl  amplitude,  as  shown  in  Figure  5,  these  equations  must  be  expressed  as 
follows: 

Eq.  3:  Yt  -  tA.COSXz,  4  Bj  SinXZi  4  C,  «sh  Xz,  4-  DySmhX^je 

Expanding  the  integration  constants  in  terms  of  their  real  and  imaginary  parts 


A, 

Thus,  T, 

where  A 

6 


A,p  ♦  i  B,  *  . 

A  cos  cert  _  B  sm  cot 

AfcosXz,  +  B,rsm XZ|  -4  C ,ra*KXz,  +  ^rs^Xz, 

4  t  *  • 

A,  COSXz,  4>  B.^S^Xz,  4  cj’  cesKXz,  4  DfsmhXZ,  131) 


Rejecting  the  imaginary  components  gives  the  above  Equations  (31) 
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Similarly,  for  the  other  coordinate  deflexions: 

Eq.  6:  X,  -  [e,  Op  Xz,  +  Fj  Xz,  4  ^,COSX  A?,  >  *'»»*  Az,  J  £ 


ilAt 


Thus 

where 


cl 

X, 

E 


£,r+  ie,1  i  £)  s  - 

E  CbS  cot  —  F  sir*  4o "t 

£,r  CoS  A2,  +■  F^SihAz,  v  ^,rCosA>2,  v  H,rsmhAz, 
E,4'  cos  A z,  +  F/smA^t  +  A  2,  + 


(32) 


The  axes  of  the  whirl  amplitude  ellipse  may  be  calculated  as  follows: 
Major  axis  : 

Minor  axis l 


a.  =  j\  [aN  tfV  eV  F‘]  +  {  J  [A\  b\  cVf*] -  4  [-  afVbe]z 

b  *  /  2  (V*  &*’♦£’>  f]- l  J  [aVbVeVf'3-H  [-AFv  be]2-  ' 


(33) 


(34) 


Angle  between  x-  axis  and  major  axis  of  ellipse: 


3 


(35) 


Angle  between  major  axis  and  unbalance  force: 


f  =  ^  rgfe;^] 


(36) 


In  the  region  0  Si  Zx  Lx  ,  rotor  displacement  amplitudes  are  given  by  Eqs, 
(4)  and  (7).  Operating  on  these  in  a  similar  manner  to  the  above  gives: 


where 


rx 

A 

6 

*z 

E 

F 


A  cos  to t  —  6  ski  cot 

ccs  AZj.  +  B^Sn  Azt  +  cj*  cosh  Xlt  +■  P^SwhAz,, 

Ai  Cos  Azx  +  bJ,  Sh  AZv  +  ci  cosh  At*  +  Sw»h  Azx 

E  cos  cot  -  F  sinust 

Ex^ATv  +  F^SnA^  V  (S^ccshA^-K  Hlsi»hA Zx 
Ex  cos  Azx  ♦  F^s/n  XT*  4  <s‘coshAzx  +  H*SmAAzx38) 


(37) 


The  major  axis  a,  minor  axis  b,  position  angle  a  and  phase  angle  j3  may  again 
be  calculated  by  using  Eqs.  (33),  (34),  (35)  and  (36). 
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Calculation  of  Transmitted  gores 


The  force  transmitted  from  the  rotor  to  the  bearings  due  to  its  motion  may  be 
calculated  by  evaluating  the  end  shear  force  acting  on  the  rotor,  as  this  force 
is  made  up  of  the  resultant  spring  and  damping  acting  at  that  point.  For  the 
left-hand  bearing  z^  =  o,  from  Eq.  (9), 


V(Y,)(o)  -  F  (Y,)(o)  »  o 


where 


-  FI  ^  3.  -££>i3|a,Si»iAz,  & 

3  L 


d.z[ 


Rejecting  the  imaginary  components,  the  dimensionless  force  in  the  y-direction 
is  given  by 


if 

w-e 

S  A  cos  cot  -  B  Sweat 

where 

A 

*  -i>  [A,rs»«Az,  -  BfcosAz,  ♦  C,r swfcAz,  +  DfcoshAz,] 

6 

*  - V  [A,1  swAz,  -  B^fccos  Az,  +  Cx  swJtAz,  ♦  CfasAAz,]  (39) 

Similarly,  the  force  in  the  x-  direction  is  given  by 


where 


5^ 

We 


E  cos  wt  -  F  SMi 


E  *  -V  [E,rs»nAz,  -  F(rto*Az,  ♦  £,rs*kAz,  ♦  HfcwAAz,] 

(40 

P  *  -y  [e,1  Sn>2,  -  cos >2,  *  ♦  h/osM*,] 

Maximum  and  minimum  values  of  the  transmitted  force  at  the  left-hand  bearing  may 
also  be  obtained  from  the  properties  of  the  transmitted  force  ellipse,  described 
by  Eqs.  (34)  and  (34)  as  follows: 


=  Ji  [av.  b*<  ev*  F “ ]  *■  i / [a*'*®*'* f‘] -  4  [-AF  »be]z 
£*♦  F*]*i  /[Abv*E1*f‘]-4[-AF+B£]* 
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For  the  right-hand  bearing  Z ^  the  transmitted  force  in  the  y-  direction  is 

found  from  Eq .  (18) 

V(XX<-0  ♦  F(Y,XO  -  o 

Thus  5l£  -  Acce,  L*t  -  B  s»w  tot 

W‘£ 

where  A  =  1?  [  >ZA  -  B^CeS  +  Cjs«nJ»Xzx  +  COSkX^  ] 

B  =  V  [  sin  Xz* -  B^coSXZx  +  4  DjfiBfcXz*]  (41) 

Similarly,  in  the  x-  direction, 

5l£  =  g  ca&ont  -  Fswtot 

W-e 

where  E  =  l)  £  sin  XZj.  -  F*d»«Xzv  •*■  >ZX  +  HcoskXZj.  ] 

F  *  V  ♦  ^SihkXz*.  *  H‘o*fcXzx]  (42) 

Maximum  and  minimum  values  of  the  transmitted  force  and  of  the  phase  angle  may 
now  be  obtained  using  the  above  constants,  by  the  previously  described  method. 

The  above  equations  have  been  programmed  for  computer  solution.  The  program  is 
given  in  Appendix  B.  Computed  results  for  maximum  amplitude  and  maximum  trans¬ 
mitted  force  are  shown  graphically  in  Figures  (6)  to  (24). 
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RESULTS 

General  Features  of  the  Results 

The  results  obtained  apply  to  ranges  of  parameters  which  identify  the  performance 
characteristics  of  a  high-speed  rotor  in  fluid-film  bearings.  Also  provided  is 
a  maximum  of  design  information.  The  location  of  the  rotor  unbalance  on  the 
shaft  determines  the  modes  which  may  be  excited  in  the  motion.  The  speed,  stiff¬ 
ness  ratio  and  bearing  eccentricity  govern  the  rotor  amplitude  and  the  transmitted 
force.  Table  II,  page  19  ,  lists  the  unbalance  location  and  amplitude  position 

or  transmitted  force  location  for  each  of  Figures  6  through  23.  Two  particular 
cases  have  been  examined  in  detail:  (a)  Midspan  (static)  unbalance  =  0.50  L 
and  (b)  unbalance  at  =  0  45  L  and  0  55  L  superimposed  180  degrees  apart  to 
give  the  dynamic  unbalance  condition.  The  charts  show: 

(a)  Variation  of  maximum  rotor  amplitude  at  a  given  station  with  speed,  due 
to  rotor  unbalance 

(b)  Variation  of  maximum  transmitted  force  between  rotor  and  bearing  with 
speed,  due  to  rotor  unbalance. 

Using  the  charts  and  principles  discussed  in  this  section,  a  preliminary  high¬ 
speed  rotor  design  may  be  established.  Final  design  behavior  may  be  examined 
by  use  of  the  computer  program  listed  in  Appendix  2.  The  following  is  a  dis¬ 
cussion  of  certain  features . 

Unbalance  at  Mid  span 


Due  to  symmetry  of  the  system,  only  symmetrical  rotor  modes  may  be  excited. 
Examples  of  these  mode  shapes  are  given  in  Figs.  26.  Amplitude  at  this 
unbalance  position  is  large  in  all  modes.  Amplitude  at  the  unbalance  position 
is  shown  in  Figs.  6  through  8,  for  eccentricity  ratios  of  0.2,  0.5  and  0.7 
respectively.  Amplitude  at  the  bearings  is  shown  in  Figs.  9  through  11.  Trans¬ 
mitted  force  at  the  bearings  is  indicated  in  Figs.  19  through  21. 

Unbalance  at  =  0.45  L  and  0.55  L,  Superimposed 

This  unbalance  condition  is  also  associated  with  large  rotor  displacements  in  the 
lowest  mode.  Although  the  unbalance  forces  are  symmetrically  arranged  about  the 
half-span  point,  the  180  degree  phase  difference  between  the  forces  causes  a 
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TABLE  II 


Details  of  System  for  Unbalance  Response  Calculations 


Fig.  No, 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


Unbalance 

Position 


Lj  =  0.5L, 
Midspan 


*  o.5L. 
Midspan 


L  =  0.45  L 


L  »  0.45L 


Amplitude 

Position 


21  "  L1 
Midspan 


Zl  =  ° 

left-hand 

Bearing 


2  0. 25Lj 


Z1  *  ° 
Left  Hand 
Bear ing 


Transmitted 
Force  Position 


Bearing 

Eccentricity  Ratio 

0.2 

0.5 

0.7 

0.2 

0.5 

0.7 

0.2 

0.5 

0.7 

0.2 

0.5 


17 


0.7 


18 

19 

20 

21 


L  *  0.45L 
L  *  0.5L 


Z1  ^  L1 


Z1  1  ° 

Left-hand 

Bearing 


0.5 

0.2 

0.5 

0.7 


22 

23 


Lj  *  0.45L 


'1 


0 


Left-hand 

Bearing 


0.2 

0.5 


24 


0.7 


moment  to  act  on  the  rotor  which  gives  rise  to  both  symmetrical  and  asymmetrical 
modes  of  oscillation.  These  modes  are  shown  in  Fig.  27.  Amplitude  response  at 
Zl  =  0.25L^  is  given  in  Figs.  12  through  14,  close  to  the  left-hand  journal  but 
sufficiently  on  the  rotor  for  the  relative  response  of  all  modes  to  be  present 
in  the  results.  Amplitude  response  at  the  left-hand  bearing  is  shown  in  Figs.  15 
through  17,  and  force  transmitted  at  the  left-hand  bearing  is  indicated  in  Figs.  22 
through  24.  The  amplitude  versus  speed  results  shown  in  Fig.  18  indicate  the 


displacement  response  at  the  unbalance  position  for  asymmetrical  unbalance 

Bearing  eccentricity  is  T}  =  0.5.  Response  for  two  rotor  stiffnesses  is  shown, 
c/5  =  0.3  and  c/5  =  10.0. 


Influence  of  Eccentricity  Ratio  t; 

Each  eccentricity  value  corresponds  to  a  particular  Sommer feld  Number  (^rp) (tb ^ 

W|  c 

The  relationship  for  the  plain  cylindrical  bearing  is  shown  in  Fig.  25.  In 
calculating  any  particular  rotor-bearing  system,  it  is  necessary  to  first  cal¬ 
culate  the  Sommer feld  Number  corresponding  to  the  bearing  operating  conditions. 
Using  Fig.  25,  the  eccentricity  can  be  found,  corresponding  to  the  bearing  L/D 
ratio.  The  eccentricity  so  found  allows  the  appropriate  figures  to  be  selected 
for  the  determination  of  amplitude  and  transmitted  force. 


Influence  of  Unbalance  on  Rotor  Dynamic  Performance 

The  rotor-bearing  system  characteristics  may  be  discussed  conveniently  in  terms 
of  the  following  seven  cases  of  response  results: 

1.  Midspan  amplitude;  midspan  unbalance.  Influence  of  eccentricity  ratio. 
Figs  6,  7,  8. 

2.  Journal  Amplitude;  midspan  unbalance.  Influence  of  eccentricity  ratio. 
Figs  9 ,  10 ,  11 

3  Shaft  amplitude  at  Z^/Lj  -  0.25;  unbalance  at  L^/L  *  0.45,  0.55. 
Influence  of  eccentricity  ratio.  Figs.  12,  13,  14. 

4  Journal  amplitude;  unbalance  at  L^/L  *  0.45,  0.55.  Influence  of  eccen¬ 
tricity  ratio.  Figs  15,  16,  17. 

5.  Shaft  amplitude  at  Z^/ =  1,  unbalance  at  L^/L  =  0.45.  Eccentricity 
ratio  T)  =  0  5.  Fig.  18. 

6.  Transmitted  force;  midspan  unbalance.  Influence  of  eccentricity  ratio. 
Figs.  19,  20,  21, 

7.  Transmitted  force;  unbalance  of  L^/L  =  0.45,  0.55.  Influence  of 
eccentricity  ratio.  Figs.  22,  23,  24. 

A  number  of  general  conclusions  can  be  drawn  from  each  case,  and  these  are  listed 
below.  The  terms  "low  speed  range,"  "medium  speed  range,"  and  "high  speed  range, 


are  used  for  convenience  and  refer  approximately  to  speed  ratios,  0  to  1.0, 

1.0  to  10.0,  and  10.0  to  24.0  respectively.  Critical  speeds  are  C 
identified  below  by  reference  to  these  speed  ranges,  because  changes  in  stiff¬ 
ness  parameter  and  unbalance  position  influence  the  modes  which  may  appear  In 
the  motion,  and  make  comparison  difficult  along  conventional  lines. 

Case  A:  Flexible  bearings  serve  to: 

(1)  increase  shaft  amplitude  in  the  low  speed  range, 

(2)  attenuate  the  fundamental  critical  amplitude, 

(3)  cause  a  critical  speed  with  large  amplitude  in  the  medium 
speed  range. 

Rigid  bearings  tend  to: 

(1)  limit  amplitude  in  the  low  speed  range, 

(2)  cause  a  critical  with  large  amplitudes  around  co/aac  »  1.0 
(Rigid  bearing  critical) 

(3)  attenuate  amplitude  in  the  medium  speed  range. 

Both  bearing  types  tend  to: 

(1)  cause  a  large  critical  in  the  high  speed  range,  with  larger 
amplitudes  in  the  case  of  flexible  bearings. 

The  effect  of  increased  journal  eccentricity  is: 

(1)  increased  amplitude  in  the  low  speed  range  for  flexible  bearing  case 

(2)  increased  the  medium  speed  amplitude  peak  for  flexible  bearings. 

(3)  increased  the  number  of  small  amplitude  peaks  throughout  the  speed 
range,  for  both  bearing  types. 

Case  B:  Flexible  bearings  tend  to: 

(1)  increase  journal  amplitude  in  the  low  speed  range, 

(2)  cause  a  large  amplitude  critical  in  the  medium  speed  range, 

Rigid  bearings  tend  to: 

(1)  limit  amplitude  in  the  low  speed  range, 
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(2)  attenuate  amplitudes  in  the  medium  speed  range. 

Both  bearing  types  tend  to: 

(1)  give  similar  amplitudes  around  a >/u>c  *  1.0, 

(2)  give  a  critical  in  the  same  area  in  the  high-speed  range,  with 
larger  amplitudes  in  the  case  of  flexible  bearings. 

The  effect  of  increased  journal  eccentricity  is: 

(1),  (2),  (3)  same  general  effect  as  in  Case  A. 

(4)  decreased  the  low  and  medium  speed  attenuation  in  the  case  of  rigid 
bearings. 

Case  C:  Flexible  bearings  tend  to: 

(1)  increase  journal  amplitudes  in  the  low  speed  range, 

(2)  give  a  large  critical  in  the  medium  speed  range, 

(3)  attenuate  amplitudes  below  a>/a>c  *  4.5. 

Rigid  bearings  tend  to: 

(1)  limit  amplitudes  in  the  low  speed  range, 

(2)  cause  a  critical  around  cu/<n  *  1.0 

c 

Both  bearing  types  tend  to: 

(1)  create  large  amplitudes  in  the  medium  to  high  speed  range, 

(2)  create  a  large  critical  in  the  highspeed  range,  with  larger 
amplitudes  for  flexible  bearings. 

The  effect  of  increased  journal  eccentricity  is: 

(1),  (3)  as  in  Case  A.  No  other  major  effect. 

Case  D:  Similar  tendencies  to  Case  B  occur. 

The  effect  of  increased  journal  eccentricity  is: 


(1),  (3)  as  in  Case  A.  No  other  major  effect. 
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Case  E:  Similar  tendencies  occur  with  bearing  transmitted  force  to  those 
which  occur  with  amplitude  in  Case  B. 

The  effect  of  increased  journal  eccentricity  is: 

(1)  Somewhat  increased  transmitted  force  in  low  speed  range  for 
flexible  bearings. 

Case  F:  Flexible  bearings  tend  to: 

(1)  increase  transmitted  force  in  the  low  speed  range, 

(2)  promote  an  amplitude  build-up  in  the  medium  speed  range, 

(3)  give  higher  over -all  transmitted  forces  throughout  the  speed  range. 

Rigid  bearings  tend  to: 

(1)  limit  transmitted  force  in  the  low-speed  range, 

(2)  attenuate  transmitted  force  in  the  medium  speed  range, 

Both  bearing  types  tend  to: 

(1)  give  similar  force  values  around  =  1.0, 

(2)  cause  a  large  increase  in  transmitted  force  in  the  high  speed 
range,  with  highest  values  for  flexible  bearings, 

(3)  sustained  high  transmitted  forces  beyond  high  speed  critical. 

The  effect  of  increased  journal  eccentricity  is: 

(1)  As  in  case  E. 

Consideration  of  the  above  list  reveals  certain  consistent  tendencies  which  appear 
in  all  cases.  These  may  be  explained  as  follows: 

(a)  increased  bearing  flexibility  tends  to  limit  rotor  amplitude  in  the  low 
speed  range.  This  is  due  to  the  normal  attenuation  which  accompanies 
any  system  flexibility  and  more  effective  use  of  squeeze-film  damping 
effects  in  the  clearance.  The  critical  speed  is  decreased  because  of 
the  added  system  flexibility.  The  upper  limit  of  the  fundamental 
critical  speed  is  the  rigid  bearing  case,  co/coc  =  1.0. 
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(b)  Attenuation  in  the  medium  speed  range  is  a  normal  condition  for  opera¬ 
tion  beyond  the  rigid  bearing  critical  speed  of  any  elastic  rotor.  In 
the  flexible  bearing  case,  there  is  some  attenuation  beyond  the  first 
low-speed  critical.  This  critical  is  predominantly  a  bearing  flexi¬ 
bility,  rigid-body  effect.  The  amplitude  peak  in  the  medium  speed 
range  is  the  second  system  critical.  It  includes  both  bearing 
flexibility  effects  and  shaft  flexure  effects. 

(c)  The  high  speed  critical  occurs  at  approximately  the  same  speed  ratio  for 
all  stiffnesses  because  here  the  bearings  are  effectively  rigid  and 

the  rotor  oscillates  as  a  pinned -pinned  beam  in  the  third  harmonic  mode. 

(d)  An  increase  in  the  number  of  small  amplitude  peaks  for  operation  at 
higher  eccentricity  ratios  occurs  because  the  respective  stiffnesses 
in  the  x  and  y  directions  of  both  bearings  become  significantly 
different  with  increase  in  journal  eccentricity.  This  tends  to  give 
rise  to  two  critical  speeds,  each  associated  with  a  particular  coordinate 
stiffness. 

Bearing  Attenuation 

Figures  19  through  24  indicate  that  greater  attenuation  of  transmitted  force  is 
achieved  by  using  a  flexible  rotor  in  relatively  stiff  bearings  throughout  the 
low  and  medium  speed  ranges. 

For  symmetrical  unbalance  the  transmitted  force  for  a  rigid-bearing  flexible- 
rotor  system  C/6  =0.3,  is  often  an  order  of  magnitude  lower  than  for  a  flexible- 
bearing  rigid  rotor  system.  For  both  systems  the  transmitted  force  is  of  com¬ 
parable  magnitude  in  the  high  speed  range,  except  at  the  high  speed  critical  peak. 
For  unsymmetrical  unbalance,  the  same  general  result  occurs,  but  to  a  somewhat 
lesser  degree. 

In  the  low  speed  range,  operation  at  a  bearing  eccentricity  of  0.5  gives  the 
greatest  attenuation  of  transmitted  force,  for  both  unbalance  conditions.  At 
this  eccentricity  the  decrease  is  in  the  order  of  15  to  20  percent  over  the  other 
eccentricities.  For  the  medium  speed  range,  the  attenuation  obtained  in  the  flex¬ 
ible  bearing  rigid-rotor  case  appears  to  increase  with  increase  in  eccentricity, 
while  in  the  rigid-bearing  flexible-rotor  case  the  attenuation  decreases  with 
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increase  in  eccentricity.  In  the  high  speed  range,  the  degree  of  attenuation  ob¬ 
tained  depends  on  both  the  nature  of  the  unbalance  and  on  the  type  of  system,  but 
the  general  result  is  that  operation  at  an  eccentricity  of  0.5  promotes  the  most 
likely  condition  for  obtaining  good  attenuation  of  transmitted  force.  In  this 
case,  specific  results  should  be  obtained  from  the  curves  themselves. 

In  general,  the  results  indicate  that  the  operating  eccentricity  at  which  the 
maximum  attenuation  of  transmitted  force  occurs  depends  on  the  stiffness  ratio  of 
the  system,  the  nature  of  the  unbalance  and  on  the  speed  of  operation.  In  the 
low  speed  range  rj  =  0.5  gives  the  greatest  attenuation,  but  for  higher  speeds 
this  simple  rule  does  not  apply,  and  the  curves  must  be  used  to  select  the  operat¬ 
ing  eccentricity,  and  to  examine  how  the  run-up  and  run-down  transmitted  force 
characteristics  will  be  influenced  by  change  in  eccentricity  with  speed. 

Comparison  With  Other  Results 

The  results  of  the  present  investigation  were  verified  against  those  given  by 
a  proven  existing  computer  program.  This  other  program  was  based  on  the .  Myklestad- 
Prohl  method  of  beam  analysis,  in  which  the  rotor  is  divided  into  a  specific 
number  of  discrete  masses,  separated  from  each  other  by  flexible  members  which 
represent  the  rotor  elasticity.  The  program  includes  the  fluid-film  bearing 
properties  as  spring  and  damping  coefficients  in  both  the  x-  and  y-directions  in 
the  same  manner  as  described  in  this  report.  The  results  of  this  comparison  are 
shown  in  Figure  28.  It  will  be  seen  that  the  agreement  is  very  close  in  all 
cases.  This  existing  program  could  have  obtained  all  the  results  given  herein, 
but  the  computing  time  involved  would  have  become  increasingly  greater  for  the 
higher  modes.  The  present  analysis  takes  the  same  computing  time  whatever  the 
mode,  and  the  results  obtained  are  exact  and  independent  of  the  number  of  discrete 
masses  into  which  the  rotor  is  divided. 

The  results  obtained  were  also  compared  with  those  given  in  the  earlier  investi¬ 
gation,  References  1  and  2.  No  satisfactory  correlation  was  obtained  in  either 
case,  for  a  number  of  reasons.  The  single  disc  rotor  of  Reference  1  has  its  mass 
and  unbalance  concentrated  at  a  single  point.  When  the  mass  is  'distributed '  across 
the  span  for  comparison  with  the  present  work,  the  unbalance  is,  in  effect,  dis¬ 
tributed  also.  This  changes  both  the  magnitude  and  nature  of  the  disturbing  force, 
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and  so  changes  the  rotor  response.  No  meaningful  comparison  was  achieved  in  this 
case.  A  similar  condition  occurred  when  a  comparison  was  attempted  with  Reference 
2.  Rotor  mass  and  unbalance  are  again  concentrated  at  discrete  points  and  so  the 
problem  of  distribution  again  arises.  This  instance  is  more  complicated  as  the 
properties  of  the  first  and  second  critical  speeds  are  obtained  uniquely  for 
either  case  from  the  equations,  by  adjusting  the  nature  of  the  unbalance  and  its 
position.  This  gives  results  which  are  properties  of  that  mode  alone.  Static 
and  dynamic  response  results  cannot  be  superimposed  into  an  overall  effect.  Both 
types  of  result  arise  naturally  in  the  present  analysis.  Also,  the  results  given 
herein  apply  for  circular  cylindrical  bearings,  whereas  the  results  of  Reference 
2  apply  to  the  150  degree  partial  arc  bearing.  In  operation  the  dynamic  proper¬ 
ties  of  both  bearing  types  are  quite  similar  due  to  the  presence  of  cavitation  in 
the  cylindrical  fluid  film.  This  means  that  the  present  resultsmay  be  used  to 
obtain  a  fair  indication  of  the  performance  of  an  elastic  rotor  in  partial  arc 
bearings  as  well. 

The  results  of  comparison  with  other  data  are  that  excellent  correlation  was  found 
where  the  rotor  and  its  operation  was  simulated  exactly,  using  a  computerized 
discrete-mass  approach  and  that  no  direct  comparison  was  obtainable  with  other 
published  work  because  an  adequate  numerical  comparison  between  the  operating 
unbalance  condition,  the  rotor  mass  distribution,  and  the  mechanism  of  modal 
excitation  could  not  be  made. 
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Mode  Shapes 
Symmetrical  Unbalance 

Figure  26  applies  to  the  operating  conditions  shown  in  Fig.  7,  C/5  =  10.0  and 
T)  =  0.5.  The  five  mode  shapes  correspond  to  major  features  of  that  curve.  Local 
amplitude  peaks  occur  at  speed  ratios  of  0.6,  2.25  and  12.0.  Local  minima  occur 
at  1.50  and  9.0.  Figure  26  indicates  che  rotor  form  under  these  conditions. 

The  damping  introduced  by  the  bearings  serves  to  limit  the  peak  amplitudes,  and 
to  increase  the  magnitude  of  the  attenuation  which  may  be  realized  between 
critical  speeds.  This  is  evident  from  the  amplitude  scales  shown,  and  also  from 
Fig.  7.  The  heavily  attenuated  low-speed  critical  at  <n/a)c  =  0.60  is  mainly 
composed  of  rigid  body  shaft  displacement  within  the  bearings,  with  a  slight 
amount  of  bending.  At  cn/^  =  1.50  the  local  maximum  attenuation  occurs,  but  the 
increased  speed  effects  plus  the  high  stiffness  and  damping  present  hold  the 
amplitude  to  two-thirds  of  the  first  crtical  amplitude  at  midspan,  and  give  rise 
to  amplitudes  at  the  ends  which  are  greater  than  those  of  the  first  critical. 

The  increased  bending  associated  with  higher  speeds  is  clearly  shown  in  the 
second  critical  speed  mode  shape  at  u)/u>c  =  2.25.  This  mode  is  similar  in  form  to 
the  fundamental  free-free  vibration  of  a  uniform  beam.  The  results  show  that 
flexible  shafts  tend  towards  pinned-pinned  beam  modes,  while  rigid  shafts  tend 
toward  free-free  modes,  because  for  a  flexible  shaft  the  bearings  control  the 
motion,  whereas  in  the  latter  case,  the  shaft  properties  alone  determine  its 
modal  form. 

A  second  attenuation  trough  occurs  at  cd/coc  =  9.00.  Amplitude  at  midspan  is 

small,  but  elsewhere  it  is  of  moderate  size  with  a  maximum  at  the  ends  of  4.2. 

The  third  critical  occurs  at  cd/cd  =  12.00.  The  motion  is  almost  wholly  bending, 

c 

and  the  form  corresponds  to  the  third  mode  of  a  free-free  beam.  The  speed- 
dependent  forcing  promotes  the  large  amplitudes  which  are  indicated. 

In  this  instance,  the  symmetry  of  the  system  and  its  excitation  do  not  allow 
motions  with  asymmetrical  modes  to  occur. 


Asymmetr ical  Unbalance 

The  influence  of  dynamic  moment  unbalance  on  the  previously  discussed  operating 


conditions,  c/5  =  10.0  and  T)  =  0.5,  is  shown  in  Fig.  27.  Figure  13  indicates  the 
amplitude  response.  The  low  speed  rigid-body  critical  speed  occurs  at  cu/co,  *  0.9. 
The  motion  is  mainly  translation  with  some  rocking  and  no  bending.  This  speed 
constitutes  a  50  percent  rise  over  the  previous  unbalance  case,  but  there  is 
significant  decrease  in  amplitude  level.  There  is  very  little  bending  present 
in  this  mode.  At  oj/coc  **  3.Q,  an  attenuation  trough  exists.  The  mode  shape  here 
is  similar  to  the  free-free  beam  fundamental.  The  second  critical  speed  occurs 
at  Wcoc  =6.0.  This  corresponds  to  the  second  free-free  beam  modal  form,  and 
the  increased  existing  force  again  begins  to  promote  large  amplitudes  with  the 
higher  speeds.  This  mode  is  almost  wholly  bending.  The  third  critical  speed 
occurs  at  a>/<nc  =  20.0,  accompanied  by  even  higher  amplitudes,  with  a  shape 
corresponding  to  the  fourth  free -free  beam  mode.  There  is  an  order  of  magnitude 
amplitude  decrease  within  the  intervening  trough  between  criticals,  which  means 
that  between  a)/a>c  =  5.0  and  25.0  midspan  amplitudes  never  decrease  below  (x/e) 

=  1.0  and  are  usually  much  larger  for  most  of  this  range.  Modal  form  within  the 
amplitude  trough  was  not  determined,  but  is  is  likely  that  it  would  be  similar  to 
the  third  free-free  beam  node,  considering  the  above-mentioned  sequence  of  free- 
free  beam  forms. 

The  moment  unbalance  tends  to  give  rise  to  even  free-free  modes  more  readily  than 
to  the  odd  modes.  These  latter  tend  to  occur  in  the  amplitude  "troughs",  which 
themselves  are  often  associated  with  large  amplitudes,  probably  because  they 
contain  suppressed  criticals. 
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CONCLUSIONS 

1.  A  study  has  been  made  of  the  dynamic  behavior  of  an  unbalanced  elastic  rotor 
supported  in  damped  fluid-film  bearings.  Extensive  numerical  results  have 
been  obtained  for  amplitude  response  and  transmitted  force,  over  wide  ranges 
of  speed,  system  stiffness  and  bearing  operating  conditions. 

2.  The  results  indicate  that  rotor  motions  are  largely  determined  by  the  inter¬ 
action  between  rotor  stiffness  and  bearing  stiffness;  and  by  the  type  of 
unbalance,  force  or  moment,  which  is  present  in  the  system. 

3.  A  flexible  rotor  tends  to  vibrate  in  modes  which  are  determined  by  the 
rigidity  of  the  bearings,  i.e.,  as  a  pinned-pinned  beam.  A  rigid  rotor 
vibrates  as  a  rigid  body  at  low  speeds,  but  where  bending  effects  predominate, 
its  motions  are  similar  to  those  of  a  free-free  beam. 

4.  Symmetrical  unbalance  excites  only  those  modes  which  are  symmetrical  about 
midspan.  Unsymmetrical  or  moment  unbalance  excites  mainly  motions  which  are 
unsymmetrical  about  midspan,  and  suppresses  the  symmetrical  motions  somewhat. 

5.  The  operating  eccentricity  corresponding  to  maximum  transmitted  force  attenua¬ 
tion  depends  on  the  stiffness  ratio  of  the  system,  the  nature  of  the  unbalance, 
and  on  the  speed  of  operation.  In  the  low-speed  range  tj  =  0.5  gives  the 
greatest  attenuation,  but  for  higher  speeds  this  simple  rule  does  not  apply 
and  the  curves  given  herein  must  be  used  to  select  the  condition  of  optimum 
operation. 

6.  The  results  obtained  may  be  used  directly  as  a  guide  towards  a  preliminary 
design.  The  performance  of  the  final  design  may  be  obtained  by  using  the 
computer  program  listed  herein. 

7.  Correlation  between  the  results  obtained  herein  and  those  given  by  a  com¬ 
parable  computer  program  of  known  accuracy  has  been  demonstrated. 
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APPENDIX  A.  Notation  for  Spring  and  Damping  Coefficients 
Reference  1  lists  the  dimensionless  bearing  stiffnesses  as: 

K  K  K  K 

xx  xv  yx  yy 

—  \o>  -  Xco  — 

c  c  c  c 

where  K  ,  K  ,  K  ,  K.  are  the  bearing  stiffness  coefficients  in  the  xx,  xy, 
xx  xy  yx  yy 

yx  and  yy  directions  respectively, 
c  is  the  bearing  clearance  and 

sw 

»  where  S  =  Sommer feld  Number 

it 

and  W,  =>  bearing  load. 

Similarly  the  dimensionless  bearing  damping  is  given  as 


cuC 

—  XjJD 
c 


_IX 

—  A JX> 
c 


cuC 

yy 


where  C  ,  C  ,  C  ,  C  are  the  bearing  velocity  damping  coefficients  in  the 
xx  xy  yx  yy 

xx,  xy,  yx  and  yy  direction  respectively. 

Notation  used  differs  bet  jen  References  1  and  2  respectively,  and  in  the  present 
report.  For  ease  of  comparison  the  following  table  has  been  prepared. 
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TABLE  II 

Comparative  Listing  of  Symbols  Used  for  Stiffness  and  Damping  Coefficients 


Lund  and  Sternlicht  (1)  Warner  and  Thoman  (2)  Present  Work 


K 

CK  /W 

K 

XX 

y 

XX 

K 

CK  /W 

K 

yy 

X 

yy 

-  K 

CD  /W 

K 

xy 

y 

xy 

-  K 

-  CD  /W 

K 

yx 

X 

yx 

C 

CB  cp/W 

C 

XX 

yy 

XX 

C 

CB  cp/W 

C 

yy 

XX 

yy 

-  c 

CB  cp/W 

C 

xy 

yx 

xy 

-  C 

CB  cp/W 

C 

yx 

xy 

yx 

APPENDIX  B.  Program  Details  And  Listing 
C  MECHANICAL  TECHNOLOGY  INC.*  LATHAM. N.Y.*  ST5-0922  J.  MICHAUD 

C  PNO 121  VIBRATION  ANALYSIS  OF  A  UNIFORM  BEAM 

DIMENSION  SPD(IOO)  .DEFL ( 1 5  )  *01  ST  ('  10  )  »  A  I  6  *6  )  *B  (  6  *6  )  *061  »D(6) 

DIMENSION  E< 12.121 *F(  12 ) .IPIVOI 12 ) .PIVOT ( 12) .ZEE1 ( 10) .ZEE21 10) 

1  READ  743 

READ  710*  SXX.SXY .SYX.SYY 
READ  710.  DXX.DXY.DYX.DYY 

2  READ  720*  NSPD *NDEFL . ND I  ST .NZEE 1 « NZEE2 » I NPUT ♦ ND I  AG . NZ 
DO  650  J J= 1 »  ND I  ST  *  5 

650  READ  730.  D I  ST  (  J J ) . D I  ST ( J J+ 1 ) . D I S T ( J J  +  2 ) *D I  ST ( J J  +  3 ) . DI ST ( JJ  +  4 ) 
DO  66 J  MM= 1 .NDEFL.5 

660  READ  730.  DEFL  (  MM ) ,DEFL ( MM+1 ) .DEFL ( MM  +  2 ) .DEFL ( MM+3  )  »DEFL (MM+4  ) 
DO  670  I  1  =  1 .NSPD.5 

670  READ  730.  SPD ( I  I ) . SPD (  I  I  + 1 ) .SPD ( I  I +2  )  »SPD ( I  I  +  3 ) . SPD ( I  I +4 ) 

IF  (NZ)  687,674.687 

674  IF  ( NZEE1 )  675,675,680 

675  READ  735.ZEE2< 1  ) 

NNN  =  2 

kFORK= 1 
ZEE=ZEE2  <  1  ) 

GO  TO  690 

680  DO  685  I=1,NZEE1,5 

685  READ  730,  ZE E 1 ( I ) .ZEE  1 < I  + 1 )  .ZEE  1  (  I +2 ) .ZEE  1 < I  +  3 ) .ZEE  1 ( I +4 ) 

IF  (NZEE2)  688.688.689 

688  NNN= 1 
KF0RK= 1 
ZEE  =ZEE  1  ( 1 ) 

GO  TO  690 

689  DO  686  I=1,NZEE2.5 

686  READ  730.ZEE2( I ) »ZEE2( I  +  1 ) .ZEE2I 1+2) »ZEE2( 1+3) .ZEE2I 1+4) 


KF0RK=2 


« 


34 


GO  TO  690 

687  READ  730.  Zl.Z2.ZEE 
KFQRK= 1 
NNN*  1 


690  WRITE 

OUTPUT 

TAPE 

3. 

762 

WRITE 

OUTPUT 

TAPE 

3. 

743 

WRITE 

OUTPUT 

TAPE 

3  . 

740. 

SXX.SXY.SYX.SYY 

WRITE 

OUTPUT 

TAPE 

3. 

741  . 

DXX.DXY.DYX.DYY 

DO  900  J  J  = 1 . ND I  ST 
DO  910  MM= 1 ♦ NDEFL 
GO  TO  (691.692) .KFORK 


691 

WRITE 

OUTPUT 

TAPE 

3. 

744 

WRITE 

OUTPUT 

TAPE 

3. 

745 

WRITE 

OUTPUT 

TAPE 

3. 

746. 

DIST  C JJ) .DEFL(MM) .NNN.ZEE 

WRITE 

OUTPUT 

TAPE 

3. 

747 

GO  TO 

695 

692 

WRITE 

OUTPUT 

TAPE 

3  ♦ 

751 

WRITE 

OUTPUT 

TAPE 

3. 

745 

695 

DO  92C 

1  11  =  1. 1 

NSPD 

GO  TO 

(810.805) .XFORK 

805 

WRITE 

OUTPUT 

TAPE 

3. 

752  ♦ 

SPD< I I ) *DIST< JJ) .DEFL(MM) 

WRITE 

OUTPUT 

TAPE 

3. 

780 

810  D I ST2= 1 .0-0 1  ST ( J J ) 

VAR  1  =  3 • 14 1 5927*SQRTF( SPD ( II)) 
ARG1=DIST( JJ)*VAR1 
ARG2=D I ST2*VAR1 

GNU=0.257021*DEFL(MM)*SPD( 1 1 )*VAR1 
VAR2*COSF( VAR1 ) 

VAR3*SINF( VAR1 ) 

VAR4*C0SF( ARG1 > 


VAR5=S I NF ( ARG1 ) 
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VAR6=C0SF (ARG2) 
VAR7=SINF(ARG2) 

VAR9  =  EXPF(VAR1  ) 

VAR8=< VAR9+1 .0/VAR9 ) / 2 • 0 
VAR9=(VAR9-1.0/VAR9)/2.0 
VAR11=EXPF IARG1 ) 

VAR10=(VAR1 1  +  1.0 /VAR  11 )/2.0 
VARll=(VARll-1.0/VARll)/2.0 
VAR13=EXPF ( ARG2 ) 

VAR12=< VAR13+1.0/VAR13)/2.0 
VAR13=( VAR13-1.0/VAR13) /2.0 
AFR5=-GNU*VAR7+SYY*VAR6 
AF  I  5=  DYY*VAR6 

AFR6=  GNU+VAR6+SYY* VAR7 
AF  I  6=  OYY*VAR7 

AFR7=-GNU*VAR13+SYY*VAR12 
AF I  7=  D Y  Y*V AR 1 2 

AFR8=-GNU*VAR12+SYY*VAR13 
AF I  8=  DY Y*VAR 1 3 

AFR9=-GNU*VAR7+SXX*VAR6 
AF I  9=  DXX*VAR6 

AFR10=GNU*VAR6+SXX*VAR7 
AFI 10=  DXX*VAR7 

AFR1 1=-GNU*VAR13+SXX*VAR12 
AF  1 1 1=  DXX*VAR  1 2 

AFR12=-GNU*VAR12+SXX*VAR13 
AFI  12s  DXX*VAR13 

THA4=VAR7+VAR13 
THA7=VAR7-VAR13 
A  ( 1  *  1 )  =2 »Q*SYY 
A  ( 1  1 2 ) *-GNU 


A ( 1 .3) =  GNU 
A ( 1 .4) =2.0*SYX 
A ( 1 »5)=0.0 
A(1 .6) =0.0 
A(2»1)=2.Q*SXY 
A (2.2) =0.0 
A ( 2 ♦ 3 ) =0.0 
A(2.4)=2.0*SXX 
A ( 2 .5 ) =-GNU 
A ( 2  * 6 ) =GNU 
A ( 3 .1 ) = VAR2-VAR8 
A ( 3 .2 ) = VAR  3 
A ( 3  «  3 ) =-VAR9 
A ( 3 .4 ) =0.0 
A ( 3  » 5 ) =0.0 
A ( 3.61=0.0 
A(4»l)=0.0 
A ( 4 . 2 ) =0.0 
A ( 4  *  3 ) =0.0 
A (4.4) =VAR2-VAR8 
A  (  4 , 5 ) =  VAR  3 
A (4.6) =-VAR9 

A ( 5.1 ) =AFR5*VAR4-AFR6*VAR5+AFR7»VAR10+AFR8*VAR1 1 
A<  5.2)=AFR5*VAR5  +  AFR6*VAR4 
A (5. 3) =AFR7*VAR11+AFR8*VAR10 
A( 5.4) =SYX*(VAR2+VAR8 ) 

A ( 5  » 5 ) =SYX*VAR3 
A  (  5.6) =SYX*VAR9 
A ( 6 . 1 ) =SXY  * ( VAR2  +  VAR8 ) 

A ( 6  ♦  2 ) =SXY  *VAR3 
A ( 6 .3 ) =SXY  *VAR9 


A (6. 4) =AFR9#VAR4-AFR10*VAR5+AFR11*VAR10+AFR12*VAR1 1 
A (6 .5 1 =AFR9*VAR5+AFR10*VAR4 

A(6»6) =AFR1 1#VAR11+AFR12*VAR10 

B  t  1  *1 ) =2.0*DYY 

B ( 1 *2 ) =0.0 

B ( 1  .3) =0.0 

B ( 1 *4) =  2  »  0*D YX 

B( 1 .5) =0.0 

B(  1  .6) =0.0 

B(2»1)=2.0*DXY 

B ( 2 .2 1 =0.0 

B ( 2  » 3 1 =0*0 

B  (  2 .4 1 =2 .0*DXX 

8(2.51=0.0 

3(2.61=0.0 

6(3.11=0.0 

B<3. 21=0.0 

B (3.31=0.0 

B(3. 41=0.0 

B ( 3.5 1 =0.0 

B ( 3.6 1 =0.0 

B<4. 11=0.0 

B(4, 21=0.0 

B(4. 31=0.0 

B (4.41 =0.0 

B (4.5 1 =0.0 

B (4,6  1  =0.0 

B  (  5.1  >=AFI5*VAR4-AFl6*VAR5t-AFI7*VAR10+AFI8*VARll 
B ( 5 .2 1 =AF I 5*VAR5+AF I6*VAR4 
B(5»3)=AF17*VAR11+AFI8*VAR10 
B  <  5.4 1 =  DYX*<  VAR2+VAR8 1 
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B  <  5,5  > =DYX*VAR3 
B(5,6)=DYX*VAR9 
B ( 6  » 1 ) =DX Y  * ( VAR2+VAR8 ) 

B  (  6 .2  ) =DXY  * VAR 3 
B ( 6  »  3 ) =  OX  Y  *VAR9 

B(6»4)=AFI9*VAR4-AFll0*VAR5+AF I  1 1  * VAR  1 0  + AF I 12*VAR1 1 
B (6.5 ) =AFI9*VAR5+AFI 10*VAR4 
B (6.6 ) =  AF I  1 1*VAR1 1+AF I  12* VAR  10 
C  (  1  )  -  0 . 0 

C  (  2  )  =0.0 
C ( 3  )  =0.0 

C ( 4 ) =0, 5*VAR1*THA4 

C ( 5  )  =0. 5*VAR1*( THA7*SYX  +IAFI6-AFI8)  ) 

C (6  )  =0 • 5*VAR 1  * ( (AFR10~AFR12)+THA7*DXY ) 

D  (  1  )  =0.0 

D  (  2  )  =0.0 

0(3)  =-0.5*VARl*THA4 
D ( 4 ) =0 • 0 

D  (  5 ) =0 • 5*VAR 1  *  < THA7*DYX  -(AFR6-AFR8)  ) 

0(6 ) =0 • 5*VAR 1  * (  (AF I10-AF I  1 2 ) - TH A 7*SX Y  ) 

IF(NDIAG)  9,15.9 
9  DO  11  1=1.6 

11  WRITE  OUTPUT  TAPE  3,  7 9 1 . A ( I »  1  )  ♦  A  (  I  , 2 ) . A ( I ♦ 3  )  .  A (  I . 4 ) » A ( I » 5 > . A ( I . 6 ) 
DO  13  1=1.6 

13  WRITE  OUTPUT  TAPE  3,  79 1  .  B  (  I . 1 ) »B C 1 , 2 ) »B < I » 3 ) »B ( I .4 ) »B (I » 5 >  ♦  B ( I  .  6 ) 
WRITE  OUTPUT  TAPE  3,  791,  C ( 1 ) »C ( 2 ) ,C ( 3 ) ,C ( 4 ) ,C ( 5 ) ♦ C < 6 > 

WRITE  OUTPUT  TAPE  3,  791.  D ( 1 ) » U ( 2 ) , 0 ( 3 > ,0 ( 4 ) .0 ( 5 ) » D l 6 ) 

C  THIb  PROGRAM  SOLVES  A  SYSTEM  OF  LINEAR  EQUATIONS  AX=B 

C  WHERE  THE  VARIABLES  ARE  COMPLEX 

f  A  =MA T R I  X  OF  REAL  PART  OF  THE  COMPLEX  SYSTEM 

C  B  =  MA T R I X  OF  IMAGINARY  PART  OF  COMPLEX  SYSTEM 


c 


C -  COLUMN  VECTOR  OF  THE  REAL  PART  OF  THE  RIGHT-HAND  VECTORS 
D=COLUMN  MATRIX  OF  THE  IMAGINARY  PART  OF  THE  RIGHT-HAND  VECTORS 


C  THE  UNKNOWN  COLUMN  MATRIX  X  IS  STORED  IN  C  AND  D 

C  C  bEING  THE  REAL  PART  AND  0  THE  IMAGINARY  PART 

15  N  =  6 

DO  10  1  =  1. N 
DO  10  J  = 1 » N 
10  E ( I » J ) = A ( I .J) 

L  L  =  N+ 1 
LM=2*N 
DO  12  1=1. N 
DO  12  J=LL.LM 
K.  =  J  — N 

12  E ( I  . J ) =-B( I »K> 

DO  14  I =LL  «  LM 
DO  14  J=1.N 
KK= I-N 

14  E(  I  »  J  )  =B  (  K.K.  •  J  ) 

DO  16  I =LL  .  LM 
DO  16  J=LL.LM 
K  =  J-N 
KK  = I-N 

16  E  (  I  ♦  J)  =AUK.K  ) 

DO  18  J= 1 » N 

18  F ( J ) =C ( J ) 

DO  22  J=LL.LM 
K  =  J— N 

22  F(J)=D(K) 

N  =  LM 

DO  20  J= 1 »N 


20  IPIV0(J)=0 
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do  550  1  =  1  .  N 
C 

C  SEARCH  FOR  PIVOT  ELEMENT 

C 

AMAX=U • 0 
DO  1U5  J= 1 » N 

IF  (IPIVO(J)-l)  60,105.60 
60  DO  100  K= 1 » N 

IF  (IPIVO(K)  -1)  80,  100.  600 
90  IF  (AMAX)  62,63.63 

62  AMAX=-AMAX 

63  IF  ( E ( J  »K )  )  64,65.65 

64  TOOT=-E < J.K  ) 

GO  TO  67 

65  TOOT=E (J.K) 

67  IF  (AMAX-TOOT)  85,100,100 
85  I ROW= J 
ICOLU  =K 
AMAX=E< J.K ) 

100  CONTINUE 
105  CONTINUE 

I  PI  VO (  I COLU )  =  I P  I  VO  <  ICOLU) +1 

C 

C  INTERCHANGE  ROWS  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 

C 

IF  ( I  ROW- 1 COLU )  140.  260,  140 

140  DO  200  L= 1 »N 
AMAX  =  E (  IROW.L  ) 

E ( I  ROW ♦ L  >  =E (  ICOLU. L ) 

200  E( ICOLU»L)=AMAX 


F4020015 
F402001 6 
F40200 1 7 


F4020030 

F4020031 

F4020032 


AMAX=F ( IROW) 


F ( I  ROW )  =F  (  ICOLU) 

F ( ICOLU ) =  AMAX 

260  PIVOT! I )=£< ICOLU. ICOLU) 

C 

C  DIVIDE  PIVOT  ROW  BY  PIVOT  ELEMENT 

C 

E ( ICOLU.  ICOLU)=1.0 
DO  350  L= 1 »N 

350  E ( ICOLU.L) =E( ICOLU.L) /PIVOT! I ) 

F ( I COLU ) *F ( ICOLU) /PIVOT ( I  ) 

C 

C  SEDUCE  NON-PIVOT  ROWS 

C 

DO  550  L 1  =  1 » N 

IF(Ll-ICOLU)  400.  550.  400 
400  AMAX=E(L1. ICOLU) 

E ( L 1 . ICOLU)  =0.0 
DO  450  L=1»N 

450  E(L1.L)=E(L1.L)-E( ICOLU.L )*AMAX 
F ( L 1 ) =F ( L 1 ) -F ( I COLU ) *AMAX 
550  CONTINUE 
600  N=N/2 

DO  641  J=1.N 
641  C ( J ) =F ( J ) 

DO  645  J  =  LL  »LM 
<  =J-N 

645  D!K)=F<J) 

CSR2=C( 1)*VAR10+C! 3)*VAR11 
CS I  2  =  0 ( 1 )*VAR10+D( 3 )*VAR11 
GSR2=C(4)*VAR10+C(6)*VARll 
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F402004B 

F4020049 

F4020050 


F4020057 

F4020058 

F4020059 


GSI2=D(4)*VAR10+D(6)*VAR11 
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DSR2  =  C< 1 )*VAR11+C< 3 )*VAR10 

DS I  2  =  0 ( 1 )*VAR11+D< 3 ) *VAR 1 0-0 . 5* VAR  1 

HSR2  =  C ( 4 ) *VAR 1 1+C ( 6 ) *VAR 1 0  +  0 . b* VAR  1 

HGI2=D<  4 )*VAR11  +  D( 6)*VAR10 

ESR2=C ( 4 )*VAR4  +  C<  5 ) *VAR5 

E5!2=D<  4)*VAR4+D(5)*VAR5 

FSR2  =  -C ( 4 ) *VAR5  +  C ( 5 ) *VAR4-0 . 5* VAR  1 

FSI2  =  -D(4)*VAR5  +  D<  5)*VAR4 

ASR  2 -C l  1 )*VAR4  +  C<2 )*VAR5 

AS  I 2  =  D ( 1)#VAR4+D(2 )*VAR5 

BSR2=-C< 1)*VAR5+C(2)*VAR4 

BSI  2  =  ~D(  1  >*VAR5  +  D<  2  )*VAR4  +  0.5>*VARl 

IF(NDIAG)  643  »644  » 643 

643  WRITE  OUTPUT  TAPE  3,  791.  C ( 1 ) »C ( 2 ) .C < 3 ) *C ( 4 J  •  C < 5  I ♦ C ( 6 ) 
WRITE  OUTPUT  TAPE  3.  791.  D (1) .0 ( 2 ) »D < 3 ) .D ( 4 ) »D ( 5 ) »D ( 6 ) 
WRITE  OUTPUT  TAPE  3.  791,  CSR2»CSI2»G5R2»GSI2  »DSR2 » DSI 2 
WRITE  OUTPUT  TAPE  3,  791,  HSR2 , H5 I  2 , ESR2 » ES I  2 , FSR2 , FS 1 2 
WRITE  OUTPUT  TAPE  3.  791,  ASR2 , AS  1 2 »BSR2 *BS  I  2 

644  GO  TO  (810,812),  KFORK 

810  IF  (NZ)  1  SOU  ,811,1 bOO 

811  GO  TO  ( 812 ,830) ,NNN 

812  DO  2000  I  =  1  »NZEE 1 
ZAR1=ARG1*ZEE1( I  ) 

VAR14=C0SF (ZAR1 ) 

VAR15=SINF (ZAR1 ) 

VAR19=EXPF(ZAR1) 

VAR 18= ( VAR  19+ 1,0 /VAR 19 1/2.0 
VAR 19= ( VAR 19- 1.0 /VAR  19 ) /2  ,0 
F I  1 =VAR 14  +  VAR 18 

AR 1 =C  < 1 )  * (  F 1 1  )+C<2)*VAR15  +  C(3)*VAR19 

B I  1 =D (  1 )  * (  F 1 1  )+D(2)*VAR15+D(3)*VAR19 


ER1=C(4)M  F  1 1  )+C(5>)*VAR15  +  C(6)*VARl9 

F 1 1 =D l 4 ) * (  F 1 1  )+D( 5 )*VAR15+D(6)*VAR19 

SUB1 1=AR1**2+BI 1**2+ER1**2+F 11**2 

SUB12=+ER1*BI 1-FI 1*AR1 

SUB13=ER1*AR1+FI1*B11 

SUB14  =  -ER1**2-FI1**2  +  AR1**2  +  BI  1**2 

SUB  15=  ER 1 *F I 1  +  AR 1  *  B I  1 

SUB16=ER1**2-FI 1**2+AR1**2-BI 1**2 

AMA XX= SORT F< SUB  1 1 **2-4. 0*SUB 12**2 ) 

AMAX1 = SORT F ( 0.5*SUB1 1+0.5 *AMAXX> 

BM I N1  =  SORT F <  0.5*SUB1 1-0.5 *AMAXX) 
0AL1=0.5*ANG(SUB14.2.0*SUB13 ) 
0BLi=0.5*ANG(SUB16.2.0*SUBl5) 

IF  (ZEEl(I)l  911.912.911 

911  AMAX=AMAX1 
BM  I N=BM INI 
OAL  =OAL 1 
OBL  =OBL  1 
GO  TO  950 

912  FAR1  =  GNU*(C(2)-C(3)  > 

FBI1=GNU*(0(2)-D(3)  ) 

FERl=GNU*(C( 5 )-C(6  )  ) 

FFIl=GNU*(D(5)-0<6)  ) 

FS  B11=FER1**2  +  FFI l**2  +  FARl**2  +  FB 1 1**2 
FS  B12  =  FER1*FBI1-FFU*FAR1 
FSB13=-FER1**2-FFI1**2+FAR1**2+FBI 1**2 
FS  B14=FER1*FAR1+FFI1*FBI1 
FMAXX* SQRTF (F SB 11**2-4.0*FSB 12**2  ) 

FMAX 1 «SQRTF (0.5*FSB11+0. 5*FMAXX ) 

FM I N1 *  SQRTF (C.5*FSB11-Q. 5*FMAXX  J 


FA1_F1  =  0.5*ANG(FSB13»2.0*FSB14) 
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AMAX*=AMAX1 
&M I N  =  BM 1 N 1 
OAL  =OAL 1 
OBL  =OBL  1 
FMAX=FMAX 1 
FM I N=FM I N1 
FALF=FALF1 
GO  TO  960 

960  GO  TO  ( 902*909) .KFORK 

902  WRITE  OUTPUT  TAPE  3.  749.  SPD( I l > .AMAX.BMlN .OAL .OBL 

GO  TO  2000 

909  WRITE  OUTPUT  TAPE  3.  76 8  . ZEE  1 < I ) . AMAX .BM I N .OAL »06L 

GO  TO  2000 

960  GO  TO  (922.928).  K.FORK 

922  WRITE  OUTPUT  TAPE  3.  770.  SPD ( 1  I > ♦ AMAX .BMI N .OAL .OBL . FMAX .FM I N. F ALF 
GO  TO  2000 

928  WRITE  OUTPUT  TAPE  3.  77 2.  ZEE  1 < I )  . AMAX . BM I N .OAL .OBL  . FMAX  .FM I N .F ALF 
2000  CONTINUE 

GO  TO  (  920.830  )  .  K.FQRK 
830  DO  3000  I = 1 .NZEE2 
ZAR2=ARG2*ZEE2(  I  ) 

VAR16  =  C0SF (ZAR2  > 

VARl 7  =  S  INF ( ZAR2  ) 

VAR21=EXPF (ZAR2 ) 

VAR20M VAR2 1  +  1.0 /VAR  21  )/2.0 
VAR21= (VAR21- 1.0/ VAR21  1/2.0 

AR2=ASR2»VAR16+BSR2*VAR17+CSR2*VAR20+DSR2*VAR21 
BI2=ASI2*VAR16+BSI2*VAR17+CSI2*VAR20+DSI2*VAR21 
ER2=ESR2*VAR16+FSR2*VAR17+GSR2«VAR20+HSR2*VAR21 
F I2=ESI 2*VAR16+FSI 2*V AR 1 7+GS I 2*VAR20+HS I 2*VAR2 1 


SUB21=AR2**2+BI2**2+ER2**2+FI2*«2 


5UB22=LR2*bI2-FI2*AR2 

SU623=ER2*AR2+F I2*BI 2 

SUb24=-ER2**2-FI2**2+AR2**2+Bl2**2 

SUb25=  ER2*FI2+AR2*bI2 

SUb26  =  ER2**2-F I 2**2+AR2**2-B  12**2 

AMAXX  =  S(JRTF(  SUB2 1 ** 2-4 . 0*SU6 2 2 * *2  ) 

AMAX2=$QRTF ( 0 . 5*5UB2 1 +0 . 5 *AMAX X ) 

BM I N2=SURTF ( C . 5*SUB2 1 -0. 5 *AMAXX ) 

OAL2=Q.5*ANG<  SUB24 . 2 . 0*SUB2 3 ) 

0BL2=0.5*ANG< SUB2 6 . 2 . Q*SUB2 5 ) 

IF  ( ZEE2 ( I )  —  1  •  0  >  914*915 *914 

914  AMAX= AMAX2 
BM I N=BM I N2 
OAL=OAL2 
OBL=OBL2 
GO  TO  970 

915  FAR2=  GNI)*<ASR2*VAR7-BSR2*VAR6+CSR2*VAR13  +  DSR2*VAR12> 
FBI  2=  GNU*( ASI2*VAR7-BGI 2 *VAR6+CS I  2* VAR  1 3  +  DS I 2*VAR 1 2 > 
FER2  =  GN0*(ESR2*VAR7-FSR2*VAR6+GSR2*VAR13+HSR2*VAR12) 
FF 12=  GNU*(ESI2*VAR7-FSI 2  *VAR6  +  GS  I  2  *  VAR  1 3  +  FiS  l  2*  VAR  1 2  ) 
FS  B21=FER2**2+FFI2**2+FAR2**2+FBI2**2 

FS  B22=FER2*FBI2-FFI2*FAR2 

FSB23=-FER2»*2-FFI2»*2+FAR2**2+FBI2**2 

FS  B24=FER2*FAR2+FFI2*FbI2 

FMAXX  =  SQRTF( FSB21»*2-4.0*FSB22»*2  > 

FMAX2=SQRTF(0.5*FSB21+0.5*FMAXX) 

FMIN2=SQRTF<  0 . 5*FSB2 1 -0. 5*FMAXX ) 
FALF2=0.5*ANG(FSB23»2.0*FSB24) 

AMAX=AMAX2 
BM I N=BM I N2 


OAL=OAL2 
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0BL=0BL2 
FMAX=FMAX2 
FMIN=FMIN2 
F  AL  F  =  F  A  L  F  2 
GO  TO  980 

970  GO  TO  ( 972*978) tKFORK 

972  WRITE  OUTPUT  TAPE  3*  749.  SPD ( I  I )  ♦  AMAX *BM I N *OAL *OBL 
GO  TO  3 UOU 

978  WRITE  OUTPUT  TAPE  3.  765.  ZEE2 (  I  )  . AMAX »6MI N .OAL .OBL 
GO  TO  300C 

980  GO  TO  (982*986).  KFORK 

982  WRITE  OUTPUT  TAPE  3.  770.  SPD < I  I >  . AMAX .BMI N .OAL »OBL . FMAX  ,FM I N . F ALF 

GO  TO  3000 

986  WRITE  OUTPUT  TAPE  3.  774.  ZEE2 ( I ) . AMAX . BM I N »OAL .OBL . FMAX »FM I N . FALF 
3000  CONTINUE 
GO  TO  920 
1500  Z AR 1 = ARG 1*Z 1 

VAR14  =  C0SF (ZAR1  ) 

VAR15  =  SINF(ZAR1  ) 

VAR19  =  EXPF(ZAR1  ) 

VAR  18= < VAR  1 9  +  1 .0/ VAR  19  )/2.0 
VAR  19= ( VAR 19- 1.0/ VAR  19  )  /2 .0 
FI1=VAR14+VAR18 

AR 1 =C ( 1  )  * (  F I  1  )+C( 2 )*VAR15+C( 3 >*VAR19 

B  1 1 =D ( 1 ) * (  F I  1  )+D( 2 )*VAR15  +  D( 3 )*VAR19 

ER1=C(4)*(  F 1 1  )-*-C(5)*VAR15  +  C(6)*VAR19 

FI1=D(4)*(  F 1 1  )+D( 5)*VAR15  +  D(6)#VAR19 

ZAR2=ARG2»Z2 
VAR16=C0SF(ZAR2) 

VAR17=SINF(ZAR2) 


VAR21=EXPF (ZAR2 ) 
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VAR20=( VAR21+1.0/VAR21 ) /2.0 
VAR21M  VAR21  — 1.0/VAR21  )  / 2  •  0 

A R2 = A 5R2* VAR 16+BSR2*VAR1 7+CSR2* VAR20+DSR2*VAR2 1 
B I2=ASI 2*VAR16  +  BSI2*VAR17  +  CSI2*VAR20  +  DS  I2*VAR21 
ER2=ESR2*VAR16+FSR2*VAR17+GSR2*VAR20+HSR2*VAR21 
f I2  =  ESI2*VAR16+FSI2*VAR17+GSI2*VAR2  0+FISI2*VAR21 
Y  Y 1 =AR 1 -AR2 


XX1-BI 1-BI2 
YY2=ER1-ER2 
XX2  =  F I  1-F  I  2 

SUB31=YY1**2+XX1*#2+YY2**2+XX2**2 
SUB32=YY2*XX1~XX2*YY1 
SUB33=YY2*YY1+XX2*XX1 
SUB34=-YY2**2-XX2**2+YY1**2+XX1**2 
SUB35=YY2*XX2+YY1*XX1 
SUB36=YY2**2-XX2**2+YY1**2-XX1**2 
AMAXX=SQRTF  <  SUB 3 1 **2-4 • 0*SU6 32**2 1 
AMAX  =  SORT F <  0.5*SUB3 1  +  0.5 *AMAXX ) 

BMIR  =SURTF(0.5*SUB31-0.5*AMAXX) 

UAL  =0.5*ANG< SUB 34 » 2 . 0+SUB3 3  ) 

OBL  =0.5*ANG(SUB36»2. 0* SUB 3 5 ) 

IF  (ZEE)  1100.1050,1100 
1050  FAR1=GNU*(C< 2 >-C < 3 ) ) 

FBI 1=GNU*(D<  2  >-D<  3)  ) 

FER 1 =  GNU* ( C ( 5 ) -C ( 6  J ) 

FFI  1 =GNU* ( D ( 5 )  -  D ( 6 ) ) 

FAR2=  GNU*(ASR2*VAR7-BSR2*VAR6+CSR2*VAR13+DSR2*VAR12> 
FBI  2=  GNU*< ASI2*VAR7-8SI 2*VAR6  +  CS 1 2 *VAR 1 3+DSI 2*VAR 1 2 ) 
FER2=  GNU*(ESR2*VAR7-FSR2*VAR6+GSR2*VAR13+HSR2*VAR12 ) 


FF 12=  GNU*(ESI2*VAR7-FSI 2*VAR6+GS I 2*VAR 1 3+HSI 2*VAR 1 2 > 


FAR=  FAR1-FAR2 


FBI  =  FBU-FBI2 


FER  =  FER1-FER2 
FF  I  =  FFI  1-FFI2 

FSB31=FER**2+FF  I  **2  +  FAR**2  +  FBI**2 
FSB32=FER*Fb I -FF I*FAR 
FSB33=-FER**2~FF  J **2+FAR**2+Fd I  **2 
FSB34=FER*FAR+FF I *FB I 
FMAXX=SQRTF ( FSB 3 1 **2-4. 0*FSB 32**2 ) 

FMAX3  =  SQRTF( 0 . 5»FSB 3 1 +0 . 5 *FMAX X  ) 

F M I N3  =  SOR  T F  (  0  •  5*FSB 3 1  -0 . 5 *FMAXX  > 

FALF  =0.5*ANG!FSB33»2.0*FSB34) 

WRITE  OUTPUT  TAPE  3 * 770 » SPD ( I  I ) »AMAX *BM I N »OAL »OBL * FMAX3 »FM I N3 *FALF 

GO  TO  920 

1100  WRITE  OUTPUT  TAPE  3*749.  SPD (III*  AMAX »  BMIN.  OAL*  OBL 
920  CONTINUE 
910  CONTINUE 
900  CONTINUE 

GO  TO  (  1  .  2  .  3500  )  .  INPUT 
3500  CALL  ENDJ06 
735  FORMAT  (E15.9) 

743  F ORMAT ( 49H  ) 

710  FORMAT ( 4E 1 5 • 9 ) 

720  FORMAT (815) 

730  FORMAT  (5E1Q.4) 

740  FORMAT! 5HGSXX  =  »E 11.5*6X4HSXY=* El 1.5.6X4HSYX  =  . El  1.5 »6X4HSYY=.E 11.5  ) 

741  FORMAT  <  5H0DXX  =  ♦ E 1 1 • 5 .6X4H0X Y  =  ♦ E 1 1 • 5  »6X4HDYX  =  *  E 1 1 . 5  *6X4HDYY  =  E 1 1.5) 

744  FORMAT! 1H1. 38  X 30HUNBALANCE  AMPLITUDE  RESPONSE) 

745  FORMAT! 1  HO  *  2 4X48HALL  VALUES  TABULATED  BELOW  ARE  DIMENSIONLESS) 

746  FORMAT! 1H0 *23X1 3HUNBAL.  Ll/L  = » F6 . 3 » 3X 1 3HST I FFN  RAT  1 0= » E 1 1 . 5  * 
12X3HZEE  *  I  2* 1H=»F6.3 ) 

747  FORMAT ( 7H0  SPEED .6X4HAMAX .8X4HBMI N ♦ 7X5HALPHA * 7X4HBETA * 9X4HFMAX . 
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18X4HFMIN  »  7X  5HPHASE ) 

751  FORMAT ( 1H1 .33X30HMODE  SHAPE  FOR  GIVEN  SPEED) 

762  FORMAT  ( lrtl ) 

752  FORMAT ( 1H021X11HSPD  RATIO  = »F6 . 2 »4X 1 1HUNBAL  L1/L=.F6.3» 

14X14HST I FFN  RATIO  =  *  E 1 1  *  5  ) 

7  80  FORMAT ( 12H02 1/Ll  22 / L2  . 4X4HAMAX . 8X4HBM I N . 5X 5HALPHA  »8X4HBETA . 8X4HFM 
1AX.8X4HFMIN.8X5HPHASE ) 

76  5  FORMAT ( 6XF  5 • 2  » 1XE1 1 • 5  » 1 XE 1 1 • 5 » 1 XE1 1 • 5  » 1 XE1 1  .5 ) 

768  FORMAT ( 1XF5.2.6XE11.5»1XE11.5.1XE11.5*1XE11.5> 

770  FORMAT! 1XF6.2.2XE 11.5. 1XE1 1.5.1 XE1 1.5.1 XE 11.5* 1XE 11.5. 1XE1 1.5.1  XEl 
11.5) 

7  72  FORMAT! 1XF 5. 2.6XE1 1.5.1 XE1 1.5.1 Xt 11.5. 1XE1 1.5.1  XE1 1.5*1  XE1 1.5.  IX El 
11.5! 

7  74  FORMAT (6XF5.2.1XE1 1.5 »1XE1 1.5. 1XE1 1.5. 1XE1 1.5 »1XE11. 5 .1XE1 1.5 .1XE1 
11.5) 

7  49  FORMAT ! 1 XF6 . 2 . 2 X E 1 1 • 5 . 1 X E 11 . 5 . 1 XE 1 1 . 5  »  1  XE 1 1 . 5  ) 

791  FORMAT! 1XE12.6»1XE12.6.1XE12.6.1XE12.6»1XE12.6*1XE12.6) 

END 
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Fig.  3  Fluid-film  Forces  on  Journal  at  Z^  «  0  and  Zg  -  Lg 


Fig.  4  Unbalance  Acting  at  »  L-^ 
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Relationships  for  Whirl 
Ellipse ' at  Unbalance 


ROTOR  AMPLITUDE  AT  MIDSPAN 
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Pig.  6  Rotor  Amplitude  vs.  Speed  -  Axially  Symmetrical  Unbalance 


ROTOR  AMPLITUDE  AT  MIDSPAN 


ROTOR  AMPLITUDE  AT  MIDSPAN 


ROTOR  AMPLITUDE  —■  AT  JOURNAL 
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Pig.  9  Rotor  Amplitude  vs.  Speed  -  Axially  Symmetrical  Unbalance 


ROTOR  AMPLITUDE  (-£-)  AT  JOURNAL 


ROTOR  AMPLITUDE  (A)  AT  JOURNAL 
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Pig.  11  Rotor  Amplitude  vs.  Speed  -  Axially  Symmetrical  Unbalance 


ROTOR  AMPLITUDE  ~  =  0.25 

6  L, 
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Mg.  18  Rotor  Amplitude  vs.  Speed  -  Axially  Asymmetrical  Unbalance 
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UNBALANCE  POSITION  ~  =0.5 
BEARING  ECCENTRICITY  if  =0.2 
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Pig.  19  Maximum  Transmitted  Force  vs.  Speeds  Axially 
Symmetrloil  Unbalance  .  . 
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Fig.  24  Maximum  Transmitted  Force  vs.  Speed  -  Axially 
' '  p  i|ffiPy4rlcal .  Onbalance 


ECCENTRICITY  RATIO  •* 
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Fig.  28  Comparison  Between  Results  of  Present  Analysis 
and  Discrete-Mass  Rotor  Calculation 
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notation 


a ,  b 


A1B1C1D1 

A2B2C2D2 


C 


c  c  c  c 

xx  xy  yx  yy 

D 


E1F1G1H1 

E2F2G2H2 


e 


e 

F 


g 

I 


i 


K  K  K  K 
xx  xy  yx  yy 

L 


V  L2 


M(X),  M(Y) 
N 

R 

S 

t 

V 

w 


Area  of  cross-section  of  rotor 
Whirl  ellipse  major  and  minor  diameter 
Constants  of  integration,  region  1 


Bearing  clearance 

Velocity  coefficients  for  bearing  damping 
Bearing  diameter 

Constants  of  integration,  region  2 

Modulus  of  elasticity 

Eccentricity  of  unbalance 

Base  of  natural  logarithms,  2.71828.... 

Force 

Gravitational  acceleration 

Second  moment  of  area  of  rotor  section 

F 

Spring  coefficients  for  bearing  stiffness 

Length  of  rotor 

Length  of  bearing 

Length  of  regions  1  and  2 

Length  ratio  L^/L 

Rotor  mass 

Coordinate  bending  moments 
Rev/Min 

Bearing  radius 

Sommerfeld  number  [“Tim2 
Time 

Shear  Force 
Rotor  weight 

W 

Static  Bearing  load,  ^ 
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w 

Weight  density 

X(i) 

Time -dependent  displacement 

x(i) 

x-displacement 

Y(i) 

Time -dependent  displacement 

y(z) 

y-displacement 

z 

Rotor  length  coordinate 

v  ••• 

Variables  defined  in  text 

a 

Phase  angle 

e 

Phase  angle 

5 

Deflexion  o|  uniformly  loaded  shaft, 

(_L)(  w£k-) 

1  384A  El  > 

n 

Bearing  eccentricity 

6 

\ 

Bearing  attitude  angle  ^  1/4 

r  w  Aco  i 

Frequency  parameter  L gif ' J 

V 

3  r  v 

Stiffness  parameter ,  (EIX.J)  ■(  — ) 

n 

3.14159  ... 

P 

Mass  density 

P 

Viscosity 

oi  Angular  velocity  rad/sec 

oi  Critical  speed  of  uniform  rotor  in  rigid 

bearings. 


